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ABSTRACT 

We present a Chandra survey of LMXBs in 24 early-type galaxies. Correcting for detection incom- 
pleteness, the X-ray luminosity function (XLF) of each galaxy is consistent with a powerlaw with 
negative logarithmic differential slope, /? ~ 2.0. However, /3 strongly correlates with incompleteness, 
indicating the XLF flattens at low-Lx. The composite XLF is well-fitted by a powerlaw with a break at 
(2.2lt° X lO^^erg s^^ and /3 =1.40^!]:}!^ and 2.84t° below and above it, respectively. The break 
is close to the Eddington limit for a I.4Mq neutron-star, but the XLF shape rules out its representing 
the division between neutron-star and black-hole systems. Although the XLFs are similar, we find 
evidence of some variation between galaxies. The high-Lx XLF slope does not correlate with age, but 
may correlate with [a/Fe]. Considering only LMXBs with Lx> 10'^'^ erg s^^, matching the LMXBs 
with globular clusters (GCs) identified in HST observations of 19 of the galaxies, we find the proba- 
bility a GC hosts an LMXB is proportional to L'^c^Fe where a = l.OI ± 0.19 and 7 = 0.33 ± O.II. 
Correcting for GC luminosity and colour effects, and detection incompleteness, we find no evidence 
that the fraction of LMXBs with Lx> 10^'' erg s~^ in GCs (40%), or the fraction of GCs hosting 
LMXBs (^6.5%) varies between galaxies. The spatial distribution of LMXBs resembles that of GCs, 
and the specific frequency of LMXBs is proportional to the GC specific luminosity, consistent with 
the hypothesis that all LMXBs form in GCs. If the LMXB lifetime is tl and the duty cycle is F^, our 
results imply ~ 1.5{tl / lO^yr)^^ F^^ LMXBs are formed Gyr~^ per GC and we place an upper limit 
of 1 active LMXB in the field per 3.4x10^ Lq of V-band luminosity. 



Subject headings: Xrays: galaxies- 
star clusters 



1. INTRODUCTION 



- galaxies: elliptical and lenticular, cD — Xrays: binaries — galaxies: 



As the end-points of stellar evolution, studies of X-ray 
binary (XRB) populations provide a valuable insight not 
only into black-hole and neutron-star demographics but 
also into the history of star-formation within ga laxies. In 
the four de cades since the discovery of ScoX-I ( Giacconi 



}pt al 



1962), generations of X-ray satellites have revealed 
a rich variety o f phenomenology in the ^300 XRBs i n 



lar interest are the X-ray luminosity functions (XLFs) of 
the XRBs, which in principle may provide vital clues as 
to the way in which the sources form and evolve (e.g. 
Belczynski et al. 2004; [vanova & Kalogera 2006). In 
star-forming galaxies, for example, the high-mass X-ray 
binary population dominates the XLF, making it much 
flatter than is typic al of LMXBs in the Milky Way or 



the Milky Way (IWhite et al.l |l995|: iLiu et al.| |2000|, bOOfl) . [Grimm et al.| |2003|) 



Unfo rtunately, etudioc of thoeo objocte are 



early - type galaxies ( Kilgard et al. 2002; Colbert et al 

Between earl y-type galax^ 



mpcrcd by ies, th e XL F appears remarkably similar (Kim & Fab 



a nuriber nf factors snrh as uncertainty in the distances 
to (and hence luminosities of) Galactic sources, the lim- 
ited numbers of sources and the mixture of old and young 
stellar populations in the Milky Way, which can make iso- 
lating different influences on the population challenging. 
Thes p problomo can, in part, bo mitigated through oboor — 



biano 2004; Gilfanov 2004), although there is some de- 
bate regarding its precise functional form. Various au- 
thors have parameterized it as a steeply falling broken 
powerlaw with a break around 2-4x10'^^ erg s~^ and a 
high-lumin osity negative logar i thmic differential s lope. 



esp 

vati ons of external , early- type galaxies , which comprise 



c± 2-3 (Barazin et alj |200C 



Kraft et al 



2001; Col 



a clean old stellar population and XRBs all at essentially 
the same distance. The absence of massive, young stars 
in these systems means that their XRB populations are 
entirely comprised of low-mass X-ray binaries (LMXBs) . 

Prior to the launch of Chandra only a small number 
of the very brightest point-sources in early type galax- 
i es could be res olved from the diffuse galactic emission 
( Fabbiano 1989). Chandra^ s, advent revolutionized this 
picture, enabling LMXBs t o be resolved in large numbers 
and studied directly (e.g. Barazin et al. _ 2001 ; iBlanton 



bert et al. 2004 ), and there may be evidence of a second 
2 X 1037 erg s-i(e.g. |Gilfanov| ^004) ). Al- 
the XLF has also been modeled as a sin- 



break at 
ternatively 

gle p o werlaw. truncated a b ove some limit ( sivakoff et al 

pool 



2003; Jordan et al. 



3t al . l boot 
2003f pabb 



t all - BOW- ■- Humphrey fe Buoto - ^M4 H ife mphroy et al 



2003|i-pabbiann.^006) for a recent rr-^rirw) Of particu- 



Kim fc Fabbiano! ( p003| ) ar- 
gued that the presence of a break in the powerlaw XLF 
of NGC 1316 can be artificially induced by source detec- 
tion incompleteness effects, after correcting for which the 
data were consistent with a single powerlaw XLF. We ob- 
tained a similar result for t he lenticular galaxy NGC 1332 
( Humphrey fc Buotc 2004). Applying an incomplet eness 
correction to a sam ple of ^^15 early-type galaxies, Kim 



^ Department of Physics and Astronomy, University of Califor- 
nia at Irvine, 4129 Frederick Reines Hall, Irvine, CA 92697-4575 



fc Fabbiano (2004) also reproduced this result, but they 
reported marginal evidence of a break at 5x 10'^^ erg s~^ 
when the data for all the galaxies were combined. 
Although the presence of the upper break is controver- 



2 
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sial, its luminosity is suggestively close to the Edding- 
ton limit of a 1.4 M© neutron star (Ledd =2-4x10'^^ 
erg s^^, depending on the composition of the accreting 



matt er and the neutron-star equation of state: Paczynsk: 



|1983D 



leading to suggestions that it may delineate the 
division between ne u tron-star and bla c k-hole accretors 
( ^arazin et al. 2001 ). Kim fc Fabbiano (2004), however, 
argued that the break in their data was at too high a lu- 
minosity to make this explanat ion tenable. Similarl y, in 



the elliptical galaxy NGC 720, Jeltema et al. (2003) re- 



ported a break at ~ 1 x 10^^ erg s~^, which they specu- 
lated may relate to a young ste llar population born in pu - 
tative recent merger activity. Bildsten & Deloye (2004) 
argued that if all bright LMXBs in early-type galaxies are 
ultra-compact binaries with He or C/O-rich donor stars, 
this might explain both the break and the low - lumin osity 
slope of the XLF found by |Kim fc Fabbiano| ( |2004| ). 

Another issue for which studies of early-type galaxies 
are useful is the role of globular clusters (GCs) in form- 
ing LMXBs. It has long been recognized that there is 



ilar to the distribution of LMXBs as a whole, and to the 
optical light. They argued this may indicate that dy- 
namical processes affecting GCs close to the centres of 
each galaxy may tr igger LMXB production. In contrast, 
Kundu et al. ( 2007| ) found that the radial distribution of 
GCs hosting an LMXB is similar to the GC population as 
a whole, and significantly flatter than the LMXB popu- 
lation, which they argued hinted at significant formation 
in the field. 

Based on small samples of a handful of galaxies, it has 
been shown that app roximately 4% of GCs are found to 
host active LMXBs (lAngehni ct alj |200l| ; [Kimdu et"ar 



2002, ^003 ; Sarazin ct al. 2003|). Using a sample of 6 



elliptical galaxies Kim et al."( 2006b| ) found evidence that 
this fraction varies considerably from galaxy to galaxy. 
There is no evidence that the properties of LMXBs in 
GCs systematically differ from those in the field. In con- 
trast, GCs which are brighter or redder are systemati- 



cally more likely to con t ain an LMXB (Angclini et al 



an o ^er-density of LMXBs per unit optical light within 
Milky Way GCs as compared to the field, indicating that 
dynamical processes within th em play an iinportant role 
in effi cient LMXB formation (Fabian et al, 1975; Clark 
|1975D . Since there are only --150 GCs in the MUky Way, 
however, there are insufficient sources to examine this re- 
lationship in detail. In contrast, m any early-type galax- 



ies ho st ric h GC populations (e.g . Gebhardt fc Kissler-, 
Pati| |1999| ; |Kundu fc Whitmore| ^00 1| ), providing ideal 



laboratories in which to inv estigate the LMXB- GC con 
nection. Using ASCA data. White et al. (2002) found a 
strong correlation between the specific frequency of glob- 
ular clusters (Sn) in a sample of 8 early-type galaxies 
and the integrated luminosity of LMXBs, inferred from 
spectral-fitting, leading them to suggest that all LMXBs 
within the galaxies were formed in GCs. With Chandra 
it has become possible to isolate individual LMXBs asso- 
ciated with GCs and to study the populations of sources 
individually. The reported fraction o f LMXBs associated 
with GCs varies from 
Humphrey 



■^20- 70% (e.g. [Angclini ct al.||2001, 
Buotd l2004| : |Kundu ct aL[|2002| ; 1x^0^ 



2005; Blanton ct al.| 2001 ). Using a sample of 4 galaxies 



sarazin ct al, (2003^ suggested this fraction varies along 



the Hubble sequence, possibly reflecting the increase in 
Sn from spiral to elliptical galaxies, and perhaps indi- 
cating different p opulat i ons o f LM XBs which fo rm in the 
field and in GCs. Irwin ( |2005| ) and |Juett| ( ^005] ) reported 



further evidence that the fraction depends on Sn- Irwin 



also argued that the total luminosity of the LMXBs is 
not strictly proportional to Sn, which would imply a sig- 
nificant fraction of sources, even in early -type g alaxies 
forms in the field. The strength of the Juett correla- 



t ion, ho wever, has been called into question by Kim et al 
( ^006b| ). 

Another way to address this question is to compare the 
spatial distrib utions of LMXBs a nd GCs, which has been 
attempted by |Kim et al| ( |2006bD and [Kundu et aL| ( |2007 , 
which was made publicly aval able after we submitted 
this present work), for small samples of galaxies (6 and 
5, respectively). However, these studies have produced 
inconsistent results, whic h may reflec t the sm all-number 
statistics of the samples. Kim et al. (2006b) found that 
the radial distribution of GCs hosting LMXBs is signifi- 
cantly steeper than that of the GC population, but sim- 



et al 1 EOOGH: Kmits ct al 



200 1|: iKundu ct al.||2002| , |2003| ; ISarazin ct al j |2003| : 



2006f|Sivakotf ct al.| p007|) 



Kim 

^ TEe 

luminosity dependence appears approximately consistent 
with the probability that a GC contains an LMXB being 
exact l y proportional to its luminos ity (e.g . parazin ct al 



p003| ; ISmits et al ||2006| ), ahhough ^ivakotf ct al.| ( |2007| ), 

based on a sample of 11 Virgo galaxies, argued that the 



imple ot ii Virgo galaxies, arg 
dependence upon mass is slightly stronger. 

The col our-dependence is i ndicative of a metallicity ef- 
fect (e.g. Kundu et al. 2003), and a number of authors 



have proposed pos sib le explanations ( f or a r eview, see 



lordan ct al. |2004| ). Maccarone ct al. ( ^004 ) proposed 
a model in which the effect arises due to the balance 
of mass transfer through irradiation-induced winds and 
Roche lobe overflow systematically varying from metal- 
poor to metal-rich systems. In this picture, LMXBs asso- 
ciated with blue GCs should exhibit sig nificant intrinsic 
absorption. However, Kim et al. ( 2006b| ) reported no ev- 
idence of spectral differences between sources in red and 
in blue GCs. Other possible explanations include the 
effect of metallic i ty on the stellar mass-radius relation 
( Bellazzini et al. |l995 , although this rt i ay no t be suffi- 
ciently large an effect: Maccarone et al. ^004 ), variation 
in the stell ar initial mass function (IM F) with the GC 
metallicity (Grindlay 1987, although see Kroupa 2002) or 
the absence ot a deep convective zone in metal-poor stars 
i mpeding mag netic braking, and hence LMXB formation 
( Ivanova 2006 , although it is not clear if this model would 
produce the observed power-law dependence of LMXB 
incidence on the GC metallicity). 

To date the LMXB-GC connection has only been in- 
vestigated with small galaxy samples, making it hard 
to draw strong conclusions about the relation between 
them in general. In particular the fraction of LMXBs 
which may be formed in the field is very uncertain. In 
addition, the presence of a break in the XLF remains con- 
troversial. The investigation of the LMXBs in a larger 
sample of galaxies is therefore vital to provide a clearer 
insight into the processes which give rise to them. In 
this paper, we investigate the LMXB population of 24 
early-type galaxies observed with Chandra and HST, fo- 
cusing on both the XLF and the LMXB-GC connection. 
This work is also part of a series investigating the X- 
ray properties of a sample of galaxies. In our previ- 
ous papers we have addressed the metal content of the 



ISM and their mass p rofiles ( Humphrey fc Buote 2006 
Hmnphrey et al. 2006). We consider here only a subset 
of the galaxies which yield interesting constraints on the 
LMXB populations. In addition, we include the nearby 
galaxy NGC 1404, the diffuse emission of which is dif- 
ficult to disentangle from that of the Fornax "cluster" 
but which has a fairly rich LMXB and GC population. 
The galaxies and a summary of their properties, and the 
Chandra observations used, are listed in Table |. In a 
companion paper, Himrphrey ( [2008 , hereafter Paper II), 
we present the HS'l' data-reduction and analysis upon 
which we build here. Of the 24 galaxies in th e sample, 
7 were included in the Kim fc Fabbianoj (2004 ) study of 
the XLF, 3 in the sample of Kim et al. f] |2006b|), who con- 



sidered the LMXB-GC connection, 4 in the smits et al 



(2006) sample, and 3 in the Barazin et al. ( 2003 ) sample. 
All errors quoted in this paper are 90% confidence limits 
unless otherwise stated. 

2. OBSERVATIONS AND DATA ANALYSIS 

For data reduction, we used the CI AO 3.3.0.1 and Hea- 
soft 6.0 software suites, in conjunction with Chandra cal- 
ibration database [Caldh] version 3.2.1. Spectral-fitting 
was conducted with Xspec 11.3.2p. In order to ensure the 
most up-to-date calibration, all data were reprocessed 
from the "level 1" events files, following the standard 
Chandra data-reduction threads^. Bad-pixel maps were 
created using the CIAO tool acisjrunJiotpix. We ap- 
plied the standard correction to take account of the time- 
dependent gain-drift and, for those galaxies observed 
with ACIS-I, we applied the standard CTI correction as 
implemented in the standard CIAO tools. 

To identify periods of enhanced background ("flar- 
ing"), which seriously degrades the signal-to-noise (S/N) 
we accumulated background lightcurves for each expo- 
sure from low surface-brightness regions of the active 
chips. We excluded obvious diffuse emission and data 
in the vicinity of any detected point-sources (see below) . 
Periods of flaring were identified by eye and excised. The 
final exposure times are listed in Table |^. 

Point source d etection was perform ed using the CIAO 
tool wavdetect (Freeman et al. 2002). Point sources were 
identified in full-resolution images of the ACIS focal- 
plane containing, where appropriate, the S3 chip and any 
other chips onto which the B-band twenty-fifth magni- 
tude (D25 ) ellipse listed in the Third Catalogue of Bright 
Galaxies ( |de Vaucouleurs et al. 1991, hereafter RC3) ex- 
tends. For AGIS-I observations, all of the ACIS-I chips 
were considered. To maximise the likelihood of identify- 
ing sources with peculiarly hard or soft spectra, images 
were created in three energy bands, 0.1-10.0 keV, 0.1- 
3.0 keV and 3.0-10.0 keV. Sources were detected sep- 
arately in each image. In order to minimize spurious 
detections at node or chip boundaries we supplied the 
detection algorithm with exposure-maps generated at en- 
ergies 1.7 keV, 1.0 keV and 7 keV respectively (although 
the precise energies chosen should make little difference 
to the results). The detection algorithm searched for 
structure over pixel-scales of 1, 2, 4, 8 and 16 pixels, 
and the detection threshold was set to 10^^, implying 
~1 spurious detection due to background fiuctuations 
per CCD. The source-lists obtained in each energy-band 

^ http:/ /cxc. harvard.edu/ciao/threads/index. html 



were combined and duplicated sources removed, and the 
final list was checked by visual inspection of the im- 
ages. For each source, we ob tained radial surface bright- 
ness profiles, as described in Humphrey fc Buotc ( 2004| ), 
which we examined to ensure that only point-like sources 
were considered (where there were sufficient photons). 
This led to the elimination of central sources in many 
of the galaxies, which were spuriously detected by the 
algorithm due to a centrally-peaked surface brightness 
profile. In a few other cases (notably in the centre of 
NGC 5846) , other "sources" eliminated in this way ap- 
peared to be associated with compact knots of X-ray 
emission, or the sharp edges of X-ray structures. The 
spiral galaxy NGC 4647 is projected onto the part of the 
image of NGC 4649, possibily contaminating our source 
list for this object. In order to minimize such contam- 
ination, we excluded all point-sources which lay within 
the D25 ellipse of NGC 4647 from further analysis. 

Spectra were extracted for each source from the 3-cr 
source regions returned by the detection algorithm. We 
note that there was very little problem of these regions 
overlapping; in the few cases where they were found 
to overlap, the regions were slightly shrunk to mini- 
mize likelihood of contamination, but this should not 
strongly affect the measured fluxes. For each source 
spectra and spectral responses were generated using the 
standard CIAO tools and local background regions were 
chosen from which background products were similarly 
extracted. The background regions were chosen to cover 
an area at least 8 times the source region, and contain- 
ing at least 50 photons, and excluding data from the 6-cr 
detection region of all detected sources. Furthermore, 
we restricted the background regions to lie entirely on 
the same chip and CCD node as the source centroid. A 
few sources were eliminated at this stage as they con- 
tained no (background-subtracted) photons; these are 
most likely spurious detections due to background fluc- 
tuations. Complete lists of the detected sources which 
lie within D25 are given in the Appendix. 

3. X-RAY SOURCE PROPERTIES 

3.1. Hardness ratios 

Individual point-sources typically contained too few 
background-subtracted photons to enable full spectral- 
fitting. Therefore to obtain a crude insight into their 
spectral pro perties, we computed spectral hardness ra - 



Kflgard et aL||2002| ; [Humphrey et al.||2002| ). 

Kl and HK2 for each source as, respectively. 



tios (see e.g. 
We defined HET 
HRl=(S-M)/(S-hM) and HR2=(M-H)/(M-hH), where S 
is the background subtracted count-rate in the band 0.3- 
1.5 keV, M is the count-rate in the band 1.5-3.0 keV and 
H corresponds to 3.0-5.0 keV. Fig |l shows, for a sub- 
sample of the galaxies, plots of the hardness ratios for 
those sources with more than 50 photons. For compar- 
ison purposes, we also show the loci in the HR2-HR1 
plane of a simple absorbed powerlaw as iVu and F are 
allowed to vary. The exact positions of these grids do 
not vary substantially between observations. 

Most of the sources are consistent with having low col- 
umn densities (for typical Galactic line-of-sight A'h val- 
ues the grid lines are practically indistinguishable from 
the A^H=0 case), and powerlaw slopes, F ^^1-3. Such 
parameters are as expected f or an LMXB populatio n in 
an early- type galaxy (see e.g. Humphrey et al. 2003| ). A 



TABLE 1 
Target list and observation log 



Name 


Typo 


Dist 
Mpc 


D25 

/ 


Lk 

1 nllT 
1O^^L0 


Age 
Gyr 


[Fe/H] 


[a/Fc] 


Obsid 


Inst 


Date 
dd/mm/yy 


Exp 
ks 


IC4296 


E- Radio-gal 


50.8 


3.8 


5.2 


12. 


-0.10 ±0.07 


0.31 ±0.05 


3394 


S 


10/09/01 


25. 


NGC720 


E5 


25.7 


4.6 


1.7 


2.9+°» 


0.17 ±0.13 


0.37 ±0.03 


492 


S 


12/10/00 


29. 


NGC1332 


S{s)0 


21.3 


4.1 


1.4 


4.1+^1 
—0.9 


-0.15 ± 0.14 


0.31 ± 0.10 


4372 


S 


19/09/02 


45. 


NGC1387 


SAB(s)0 


18.9 


3.3 


0.78 






4168 


I 


20/05/03 


45. 


NGC1399 


cD;Elpec 


18.5 


6.9 


2.1 


12. ± 2. 


-0.11 ± 0.08 


0.37 ±0.03 


319 


S 


18/01/00 


56. 


NGC1404 


El 


19.5 


4.1 


1.5 


12.+2- 

— 3. 


-0.28 ± 0.09 


0.25 ± 0.03 


2942 


s 


13/02/03 


29. 


NGC1407 


EO 


26.8 


5.3 


3.1 


12. ± 1. 


-0.18 ± 0.05 


0.33 ± 0.01 


791 


s 


16/08/00 


40. 


NGC1549 


EO-1 


18.3 


4.7 


1.3 


5.1 ± 0.4 


—0.10 ± 0.04 


0.240 ± 0.006 


2077 


s 


08/11/00 


22. 


NGCl 553 


SA(rl)0:LINER 


17.2 


5.3 


1.9 


5 7+°-'' 

■^•'-0.5 


23 + 02 


1 70 + 006 


783 


s 


02/01 /OO 


14. 


NGC31 15 


SO 


9.00 


7.3 


0.74 


15. ± 3. 


030 + 10 


1 1 + 07 


2040 


s 


14/06/01 


36. 


V_X ij !_> O fJ 


E7/S0 


18.6 


6.1 


1.5 


10 +3- 


_n 1 Q -1- f) 1 1 


n on + n 07 


2078 


g 


03/06/01 


35. 


NGC3607 


SA(s)0 


21.2 


4.5 


1.5 


15. ± 2. 


—0.13 ± 0.11 


0.19 ± 0.07 


2073 


I 


12/06/01 


38. 












q+0.4 

0.0_Q 2 


n 1 -t- n f\f. 


n "^A 4- n 09 


lODO 


Q 



1 A /OR /ni 


0.0 


NGC4125 


E6p6C5Liner 


22.2 


5.8 


1.8 


13. ± 5. 


—0.39 ± 0.17 


0.33 ± 0.10 


2071 


s 


09/09/01 


63. 


NGC4261 


E2-3;Liner;Sy3 


29.3 


4.1 


2.2 


15.0 ±0.6 


-0.21 ±0.06 


0.25 ±0.01 


834 


s 


06/05/00 


34. 


NGC4365 


E3 


19.0 


5.8 


1.6 


q+5-9 
j.y_Q 7 


-0.020 ±0.12 


0.19 ±0.09 


2015 


s 


02/06/01 


40. 


NGC4472 


E2/S0(2);Sy2 


15.1 


9.7 


3.2 


9.0 ± 1.2 


-0.0100 ±0.12 


0.16 ±0.02 


321 


s 


12/06/00 


32. 


NGC4494 


E1-2;LINER 


15.8 


4.5 


0.81 


15. ± 5. 


-0.24 ±0.12 


0.16 ±0.08 


2079 


s 


05/08/01 


15. 


NGC4552 


E;LINER-HII 


14.3 


5.0 


0.85 


12. ± 1. 


-0.050 ±0.061 


0.24 ±0.01 


2072 


s 


22/04/01 


54. 


NGC4621 


E5 


17.0 


5.0 


1.2 


7 -,+4.8 
'■^-2.3 


-0.030 ±0.11 


0.28 ±0.09 


2068 


s 


01/08/01 


25. 


NGC4649 


E2 


15.6 


7.4 


2.4 


13. ± 1. 


0.050 ± 0.085 


0.25 ±0.01 


785 


s 


20/04/00 


21. 


NGC5018 


E3 


42.6 


3.6 


3.0 




0.15 ±0.17 


0.0100 ±0.085 


2070 


s 


14/04/01 


28. 


NGC5845 


E 


24.0 


0.90 


0.27 


12. 


-0.12 ±0.16 


0.26 ±0.10 


4009 


s 


03/01/03 


30. 


NGC5846 


E0-1;LINER-HII 


21.1 


3.8 


1.5 


15. ± 1. 


-0.18 ±0.10 


0.22 ±0.01 


788 


s 


24/05/00 


23. 



Note. 



zero-point (Jensen et al 



Targpt lisf. a 



piijibservation details. 



All distances (pist^ arp takori froTn Tonry et al. (2001), corrected for the ttip npw r!j»plipjrl 
(pOOa), and NGG50lVffom Dn-CT relation ( [Faber et a\\ 



2003), except for IC4296, taken from Jensen et al, 
and twenty-htth magnituae isopliote (D25) diameters are taken trom Klj3' 
(20001), adopting Mkq = 3.41. We report the age of the stellar population, its mctallicity, [Fe/H] and its «-to-Fe ratio, [«/Fe] (sec [6.1 



K-band luminosities (Lk) were taken from UMA 
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their l-CT errors. The Chandra observation identifier is given (ObsID), as is the ACTS instrument (inst), the date of the observation start and 
the total exposure (exp), having cleaned the data to remove flaring (see text). 



small number of sources exhibit much harder colours, im- 
plying intrinsic absorption. Some or all of these objects 
may be heavily-absorbed background AGN. We do not 
exclude such objects from subsequent analysis, since we 
explicitly take account of the expected numbers of back- 
ground sources. 

3.2. Fluxes 

Since direct spectral fitting was impractical, it was nec- 
essary to adopt a canonical model for the spectrum to 
obtain a flux estimate for each source. We adopted a 
simple powerlaw model, with T = 1.55, which has been 
shown to fit adequately the composite spectra of detected 



estimating the flux, and so the results are self-consistent. 
Given the narrow band in which we compute fluxes, by 
definition, they underestimate the bolometric flux of the 
source. However, the extrapolation of the spectrum out- 
side the adopted band is highly uncertain and so we com- 
pute the XLF of the sources for this narrow band. The 
luminosities are shown for each source in the Appendix. 

3.3. Composite spectra 



snown to nt aaequatel y 
LMXB ( llrwin et al^ p 



2003 



sec also § 3_^), and modi- 
fied by photoelectric ab sorption due to the Galac tic ISM 
along the line-of-sight (Dickey & Lockman 1990). Such 
a model is b roadly consistent with the source hardness 
ratios (§ ^]^) and models of a similar shape have also 
been fitted to the spectra of LMXBs within th e Milky 
Way (e.g. [Church fc Balucihska-Church| p001| ). This 
model was folded through the response matrices gener- 
ated at each source position, and used to infer a counts- 
to-flux conversion ratio in the 0.3-7.0 keV band. The 
0.3-7.0 keV background-subtracted fluxes of each source 
were extracted from the regions defined from which to 
extract spectra. Although a s mall number of sources ap- 
pear heavily absorbed (§ 3.1), and so this count-to-flux 
conversion ratio is not formally correct, the X-ray lumi- 
nosity function we adopt to account for interlopers (§ [l.2| ) 
similarly does not incorporate the intrinsic absorption in 



With a sample of 15 galaxies, Irwin et al. ( pOOSj) showed 
that the composite spectrum of all the detected LMXBs 
is featureless and can be well-fitted by simple empiri- 
cal models such as a powerlaw with F =1.56 or ther- 
mal bremsstrahlung with kT=7.3 keV. In order to test 
this with our sample, for each galaxy, we accumulated 
composite source and background spectra for all of the 
observed sources which lie within D25 , excluding the cen- 
tral ~20" (where background subtraction may be more 
problematical). These spectra were added in count-rate 
space, having first scaled the background spectra by 
the ratio of the source to background BACKSCAL key- 
words. The response matrices were averaged together 
with the Heasoft tools addrmf and addarf and adopt- 
ing the background-subtracted count-rate of each corre- 
sponding spectrum as a weighting factor. 

We fitted the composite spectrum in the range 0.5- 
7.0 keV, having first rebinned it to ensure S/N of 3 in 
each bin, and at least 20 photons (to allow the use of the 
fitting statistic). We fitted the spectrum using Xspec. 
The spectrum is remarkably featureless (Fig |^) and well- 
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Fig. 1. — Hardness ratio plots of point sources with more than 50 photons for a subsample of the galaxies. The definitions of the hardness 
ratios HRl and HR2 are given in the text. Also shown is a grid of loci in the colour-colour plane for a simple photo-absorbed powerlaw as 



A'^H S'lid r are allowed to vary. TVjj increases to the left and grid lines are shown for Ah = 
increases upwards, with grid lines shown for F =0.0-3.0, in steps of 0.5 



:0.0, 0.25, 0.5, 1.0, 2.0, 4.0 and S.OxlO^^cm" 
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TABLE 2 

LMXB SPATIAL DISTRIBUTION GOODNESS-OF-FIT 



Energy (keV) 

Fig. 2. — Composite spectra of the LMXBs (upper spectrum; 
black), folded through the instrumental response. The data are 
shown fitted with an absorbed bremsstrahlung model, with kT 
fixed at 7.3 keV. Also shown are the composite spectrum of LMXBs 
associated with red GCs and, scaled for claiity by a factor 0.3, that 
of LMXBs associated with blue GCs (§5.4.2). These spectra are 
shown fitted with the same bremsstrahlung model. 



Galaxy 


prob (light) 


prob(GC,tot) 


prob(GC, weight) 


NGC720 


0.35 






NGC1332 


0.10 


0.15 


0.090 


NGC1399 


1.0 X 10-^ 


0.48 


0.85 


NGC1407 


0.20 






NGC1549 


0.12 






NGC3115 


0.78 


0.010 


0.060 


NGC3923 


0.20 






NGC4125 


0.14 






NGC4365 


0.070 


0.070 


0.22 


NGC4472 


0.23 


0.29 


0.97 


NGC4552 


0.080 


0.010 


0.44 


NGC4621 


0.24 


0.13 


0.21 


NGC4649 


3.0 X 10-3 


0.32 


0.44 



Note. — The probability that the radial distribution of 
LMXBs is the same of the optical light (prob(light)) in each 
galaxy and the probability that it is the same as that of the 
GC population (prob(GC,tot)), or a weighted combination of 
the red and blue GC distributions (prob(GC, weight); see text). 

is in better agreement with the average Galactic line- 
of-sight column densit y for the galaxies in our sample 



fitted by a single, absorbed powerlaw (x^/dof= 



with r = 1.68±0.04 and iVH=(7.5 ± 1.1) x lO^Ocm"^. If 



354/339), 

2 



we fitted the data, instead, with a bremsstrahlung model, 
we obtained a similarly good fit (x^/dof=357/339), with 
kT=7.8t°-? and A/'h=(2.1 ± 0.08. 

( |2003| ) 



10^°cm ^. These re- 



sults are consistent with [rwin et al 



although r 

is slightly larger for the powerlaw case. It is interesting 
to note that A'h obtained for the bremsstrahlung case 



(~ 3 X lO'^^cm"^) (see [Humphrey fc Buote|^006D , imply- 
ing it is more representative. I'itting this model to the 
composite source spectra accumulated for each galaxy 
individually, with Nn fixed to the nominal value for the 
appropriate Galactic line-of-sight, we obtained accept- 
able fits for all of the galaxies. 

3.4. Spatial distribution 
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Fig. 3. — Spatial distribution of the LMXBs in those galaxies with sufficient sources. Shown are the data-points having corrected for 
spatial variation in source detection incompleteness (circles; black) and, where appropriate, the best-fitting model to a suitably weighted 
combination of the distributions of red and blue GCs (solid line; red). The distribution of the GCs, corrected for incompleteness, is also 
shown, with arbitrary scaling for clarity (crosses; orange) and the same best-fitting model. The dashed line (blue) shows the distribution 
of the optical light. For comparison, we also show the data-points having not corrected for the spatial variation in source incompleteness 
(triangles; light blue). 



We next investigated the spatial distribution of the 
point-sources in a subset of the galaxies which contained 
sufficient sources. We created a set of radial bins centred 
on the X-ray centroid and each containing at least 15 
sources. We compared the radial profiles of the sources 
with the optical light, which we modeled as a de Vau- 
couleurs profile, the effective radius of which was fixed to 
t he appropria te K-band value determined from 2MASS 
(Jarrett 200Cl| ). We also included an additional, con- 
stant, component to account for interloper sources, which 
should be approximately uniformly sprinkled over the 
Chandra field. Fitting was performed with a Cash-C 
statistic mimimization algorithm. 

Since the effects of source detection incompleteness de- 
pend on the density of sources, the width of the point- 
spread function and t he "background" count-rate (e.g. 
Kim fc Fabbiano| 2003), it varies spatially across the field 



of view. It was therefore necessary to take this into ac- 
count when fitting the radial profile of sources. We dis- 
cuss in general how we estimated the effects of source 
detection incompleteness on the XLF in § 4.1. We com- 
puted a similar correction appropriate for each radial bin, 
with which we weighted an assumed X-ray luminosity 
function (assumed t o t he best-fitting broken powerlaw 
XLF we found in § |4.2D that was then integrated from 
10^^-2 X lO^^erg s~ so as to estimate the fraction of 
sources detected as a function of radius. 

In Fig H we show the spatial distribution of the point- 
sources in the subsample, both with a correction for the 
spatial variation in completeness applied, and with no 
correction. Outside the innermost 0.5' (in which incom- 
pleteness effects are most severe) the radial distribution 
of sources agrees very well with the optical light. Inside 
this region we find an excellent agreement between the 
source distribution and the optical light in approximately 
half of the tested galaxies. In the remaining galaxies, 
the optical light appears to over-predict the number of 
sources considerably. To determine the goodness-of-fit, 
we performed 100 Monte-Carlo simulations (increasing 
the number if the estimated goodness-of-fit was below 
1%), where artificial data-sets were simulated from the 
best-fitting models and were fitted with the same model. 
The fraction of simulations with a larger Cash-C statis- 
tic (i.e. poorer fit) than for the real data was adopted as 
the null hypothesis probability (prob(light)), as shown in 
Table |. If the LMXBs are, in general, distributed like 
the optical light, the distribution of prob(light) should, 
for these 13 galaxies, be uniform over the interval 0-1. 
U sing a Kolmogorq v-Smirnov (K-S) test (as implemented 



. Press et al. 1992), the probability the distributions are 
the same is only 2 x 10~^. Considering only those galax- 
ies for which we also have sufficient-quality GC data to 
measure the GC radial distribution as well (below), this 
probability is 0.05%. 

If a substantial fraction of LMXBs are produced in 
GCs, we might expect the radial distributions of the 
two populations to be similar. Since the field of view of 
the WFPC2 is significantly more restricted than Chan- 
dra, it was necessary to fit the LMXB and GC data 
with the same model simultaneously. To parameter- 
ize the GC distribution, we adopted a "beta" model 
(oc (1 -I- {R/ RcYY''^^^'^ , where Rc and e are parameters 
of the fit), and we additionally allowed a constant back- 
ground component for the LMXBs (as discussed above). 



The data were of sufficient quality to enable this for 8 
of the galaxies, and the null hypothesisorobabilities are 
reported for this comparison in Table We find mod- 
erately good agreement in most of the galaxies between 
the two distributions, although in two galaxies the null 
hypothesis probability was <5%. Blue GCs are known 
to have a considerably flatter distribution than their red 
counterparts (e.g. Paper II), while LMXBs are preferen- 
tially found in red GCs, and so a comparison with the 
radial distribution of the total GC population may not 
be entirely appropriate. We therefore also compared the 
radial distribution of the LMXBs to a weighted combi- 
nation of the radial profiles of the red and blue GCs, 
taking into account the fraction of LMXBs expected to 
be produced in each population. For this model we found 
reasonably good agreement between the distributions in 
all of the galaxies. The null hypothesis probabilities for 
these fits are shown in Table S. Based on a K-S test, 
as described above, we found that the distributions of 
these probabilities agrees with expectation (a uniform 
distribution) with a null hypothesis probability of 17%, 
significantly higher than the equivalent comparison with 
the optical light (0.05%). 

4. X-RAY LUMINOSITY FUNCTIONS 

4.1. Source detection incompleteness 

In order to measure the point-source XLF it is ex- 
tremely important to take into account the effects of 
source detection incompleteness and the Eddington bias, 
wher e the flux measurement e rr ors distort the XLF shape 



(e.g. |Kim fc Fabbiano| |2004D . |Kim fc Fabbiano| ( |2003| ) 
presented a possible strategy to correct the data for these 
effects. In order to maintain the statistical integrity of 
the data, we preferred to apply a correction to the model, 
and so adopted a modified strat egv, which in corporates 
some aspects of the algorithm of Wang] ( [2004 ) . 

For each galaxy, we performed a set of Monte-Carlo 
simulations, initially using 20 flux bins (corresponding 
approximately to equally-spaced bins from 1-20 detected 
counts), which entailed adding simulated sources to the 
image of each galaxy and then assessing the performance 
of the detection algorithm at finding and characterizing 
them. For many of the galaxies, we found that almost all 
of the sources expected to contain ^^20 photons are de- 
tected, and so we assumed that any brighter sources will 
always be detected. In those galaxies for which this was 
not the case, we added additional simulations, logarith- 
mically spaced in flux space, until ~95% completeness 
was achieved in any bin. In each flux bin, we performed 
20 such Monte-Carlo simulations, in each of which we 
added 20 sources; adding them 20 at a time does not sig- 
nificantly increase source confusion since most galaxies 
had many more than 20 LMXBs in the total field of view. 
We assumed that the point-sources were distributed ap- 
proximately as the optical light (§ 3^). For each sim- 
ulated source, to determine the number of photons to 
simulate we computed a count-to-flux conversion by cre- 
ating local ancillary spectral response files and an on-axis 
redistribution matrix and folding our preferred powerlaw 
spectral model through them. The point-source images 
were derived at each point by bilinear interpolation of 
the Caldh PSF images for a 1 kcV point-source at var- 
ious focal-plane positions, and finally Poisson noise was 
added. 



TABLE 3 

XLF FITTING RESULTS 



Galaxy N^rc f cover Lx,min Lx,20 Prob beta Lx prob(Ho) /3l Lbreak Lx 
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Note. — The best-fitting parameters derived fron 
show fit results for the composite XLF of LMXBs i: 
XLF (Nsrc), the fraction of the total galaxy light v 
contributing to the XLF (L x min) ' ^^"^ approxin- 



fitting the XLF of each galaxy, and the composite XLFs of all the sources ("composite"). We also 
L GCs ( "GC" ) and those in the field ( "field" ) . Listed are the number of sources used in fitting the 
ithin the region in which the XLF is computed (f cover )i the approximate minimum L^ of sources 
ate Lj^ of a source containing 20 photons (,L 20^' the null hypothesis probability (prob(HQ )), 
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In any arbitrary region of the Chandra field of view, 
for each flux bin this yielded the fraction of simulated 
sources which were detected and the measured count- 
rates (hence fluxes) of these sources. For brighter sources 
than those simulated, we assumed 100% completeness 
and that the distribution of measured fluxes is given by 
the Poisson statistics. Since the XLF comprises a his- 
togram of the numbers of sources detected in a series 
of narrow flux bins, this information is sufflcient to con- 
struct a "response matrix" through which any model to 
fit the XLF can be fold ed, analogous to X-ray spectral- 
fitting (e.g. Davis 200l| ) to take account of both incom- 
pleteness and the Eddington bias. We note that this 
method also applies the same correction to the interlop- 
ers, which is formally incorrect since their spatial distri- 
bution is much flatter than the optical light. However, 
there were few interlopers in the regions of interest so 
this should not significantly affect our results. 

4.2. XLF fitting 

To measure the XLFs of the galaxies, we computed his- 
tograms of the measured numbers of sources as a function 
of flux. We considered only point-sources within D25 and 
excluded the innermost 20" (where source detection may 
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be uncertain; § 3.4). The data were regrouped to ensure 



at least 5 sources per bin, and the differential XLFs fltted 
using dedicated software which folded the model through 
the response described in § 4.1. To perform the flt we 



adopted the Cash-C statistic and parameter space was 
searched for a global minimum. To assess the goodness- 
of-flt, we adopted a Monte-Carlo strategy. Adopting the 
best-fitting model, we simulated 100 fake data-sets which 
were subsequently fltted. We adopted the fraction of sim- 
ulations which yielded a best-fltting Cash statistic value 
higher (i.e. worse) than obtained for the real data as the 
null hypothesis (the model describes the data) probabil- 
ity (prob(Ho)). To account for interlopers, we adopted 
the hard-band XLF relation for sources in the Chandra 
Deep Field- South (|Tozzi ct alj^00l|) , correcting to the 
spectral-fltting band used in the present work, and the 
area of sky under scrutiny. Since the slope of the inter- 
loper XLF at high lum inosities is similar to th at of a typ- 
ical early- type galaxy ( Kim fc Fabbianc 2004 ) , it was not 
possible to flt the background normalization freely and 
so it was fixed. To investigate the sensitivity of our re- 
sults to this assumption, we have experimented with the 
effect of varying the normalization by ±50%. We found 
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Fig. 4. — Individual differential XLFs lor half of the galaxies (continued in Fig shown with the best-fitting power law model (solid line; 
red), corrected for source detection incompleteness and the Eddington bias. The expected contribution from background AGN is shown as 
a dashed line (blue-green) as is the best-fit model to the composite XLF of all the galaxies (dash-dot-dot-dot; blue). For comparison, the 
best-fit simple power law model, without incompleteness correction, is also shown (dotted line; magenta). 



that this typically caused changes in the best-fitting pa- 
rameters which were smaller than the statistical errors. 

We found that we were generally able to fit the XLF 
of each galaxy adequately with a simple powerlaw model 
of the form: 



dN 

11 



(X L 



-/9 



(1) 



The best-fitting results are shown in Table ||, and the 
XLFs arc shown in Figs M. In Table | we list the total 
luminosity in the range 10 -2 x lO^^erg s~^ determined 
from our fit, having been corrected upwards by the frac- 
tion of optical light which falls outside the region in which 
the XLF was computed. This correction allowed us ap- 
proximately to take account of the fact that not all the 
LMXBs are expected to lie in this region, since the spa- 
tial distribution of the LMXBs is generally close to that 



of the optical light. For NGC 5845 there were insufficient 
sources within the fitting region, and so we omit it from 
this table. For the 7 galaxies in o ur sa mple which over- 
lapped that of Kim fc Fabbianci| (2004), we found good 
agreement, within th e error -bars, between the measured 
values of [3. Xu et al. ( 2005 ) obtained an incompleteness- 
corrected fit for NGC 4552, but found a slightly steeper (3 
(~2.2) than we obtained. The reason for the discrepancy 
is unclear, although these authors fitted the cumulative 
luminosity function which is difficult to interpret due to 
correlations between adjacent data-points. Given that 
a single powerlaw model was typically an adequate fit, 
we found that we were generally unable to constrain the 
parameters of a broken powerlaw fit to the data. 

There was remarkably little variation in /3, indicating 
striking uniformity in the XLF shape of LMXBs in early- 
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Fig. 5. — Same as Fig M for the remaining galaxies. 



type galaxies. However, there is evidence of a small, but 
statistically-significant variation in slope. In Fig ^ we 
show how f3 varies as a function of the X-ray source detec- 
tion completeness, f^, defined as the fraction of actual 
sources with Lx= 10^^-2 x 10'^^ erg s^^ which are de- 
tected. This is obtained by integrating the completeness- 
corrected and the uncorrected XLFs over the appropriate 
flux range. To prevent spurious correlations, we fixed 
f3 — 2.0 for this calculation. The data show a clear 
anti-correlation; using Spearman's rank-order correlation 
test, we found a probability of ~ 7 x 10"''^ that the 
data are uncorrelated. This correlation indicates that the 
XLF, on average, must flatten at low Lx- If the "true" 
XLF is, in fact, a broken powerlaw which flattens be- 
low the break then the measured /3 will be an "average" 
of the slopes above and below it. As the data become 
more incomplete there are fewer sources below the break 
to weight the fit, so /3 asymptotes to the high-Lx slope. 



Since /3, where the data are very complete, appears to 
asymptote to a value significantly fiatter than 2.0, we 
conclude that /3 at low luminosities must, similarly, be 
significantly less than 2.0. We note that, if we exclude 
the two most complete data-points from this compari- 
son, we find that the probability of no correlation is still 
~ 2 X 10^^, indicating that this result is not solely driven 
by these data. 

4.3. Composite XLF 

Given the trend for systematically flatter XLFs with 
increasing source detection completeness, we next exper- 
imented with adding the XLFs in order to improve S/N. 
We combined all detected sources for all of the galaxies 
into a single source list, and used that to generate the lu- 
minosity function data-points. To account for source de- 
tection incompleteness and the Eddington bias, we per- 
formed a weighted addition of the individual responses 




Fig. 6. — Plot of XLF slope, /3 as a function of source detection 
completeness, f^. Note the correlation. 




1 10 
Lx (10^° erg s"') 

Fig. 7.— Composite XLF of all LMXBs, fitted with a simple 
power law (dashed line) and a broken power law (solid line). The 
models have been corrected for source detection incompleteness 
and the Eddington bias. 

generated for each galaxy, using the number of observed 
sources as the weighting factors. We grouped the data 



4.2, and similarly adopted the Cash-C 



identically to 
fit-statistic. 

We first fitted a single power law to the data, includ- 
ing a component to account for the expected background 
sources (estimated by adding the contributions antici- 
pated for each galaxy). We found this model to be a 
very poor fit to the data (Fig The best-fitting /3 (1.7) 
was much lower than we obtained for each galaxy in- 
dividually. This is easily understood since this shallow 
slope is largely driven by the apparent paucity of sources 
at low Lx (< 6xl0'^''erg s~^), consistent with the XLF 
fiattening at low Lx- We next tri ed to fit the da t a wit h 
the truncated powerlaw model of ^ivakoff et al. (2003), 
which gave a break at 6.7xl0'^®erg s^^ but failed to re- 
produce very well the shape of the XLF at low lumi- 
nosities (prob(Ho)< 0.1%). Finally, we fitted a broken 
powerlaw to the composite XLF. This gave a significantly 
better fit to the data; although the formal prob(Ho) is 
only 3%, the deviations from a good fit are primarily due 
to individual discrepant data-points which are unlikely 
to bias the fit. The best-fitting parameters are shown in 



Table ^, and the best-fit model is shown in Fig ^ 

We caution that, at the lowest luminosities the 
composite XLF is largely dominated by two galaxies, 
NGC3115 and NGC4552 and there exists a possibility 
that the XLFs of these two systems are not representa- 
tive of the whole population of early- type galaxies. If 
that is the case, it may bias the shape we find for the 
composite low-Lx XLF. With our current data it is not 
possible entirely to rule out such a possibility; instead it 
must be tested using very deep Chandra observations of 
other galaxies. In the present work we take this fit at 
face value, but are careful to investigate the sensitivity 
of any derived quantities to the low-Lx slope. 

To ensure that the broken power law fit did not arise 
as an artifact of our analysis, we experimented with sim- 
ulating a simple power law XLF appropriate for each 
galaxy, with f3 = 2.0. Combining these identically to the 
real data, we found that a simple power law model gave a 
good fit to the data. Conversely, if we simulated simple 
power law XLFs but with f3 fixed to the best-fit value 
for each galaxy, the composite XLF was very poorly fit- 
ted by a simple power law, and instead a broken power 
law was required. However the implied break luminosity 
(~ 1.3 X lO'^^erg s~^) was considerably higher than that 
of the composite XLF and the best-fit model for the real 
data was an extremely poor fit. 

Since the response matrices we used to correct for the 
effects of incompleteness were generated by adding a fi- 
nite number of artificial sources to each image, this intro- 
duces a statistical error into the matrix, which we have 
not explicitly taken into account in our analysis. In or- 
der to assess whether the magnitude of this effect can 
lead to significant errors in our measured XLF shape, we 
performed 100 Monte-Carlo simulations, in which all the 
"response matrices" were re-created, then combined and 
used to fit the existing data. In the interests of speed, 
we did not re-create the matrices from first principles (as 
described in § 4.1), but instead, for each luminosity bin 
of artificial sources which were added to the images, we 
drew 100 random numbers, which were distributed in the 
same manner as the observed luminosities in this bin (in- 
cluding the fraction of non-detections). Treating these 
random values as though they were the recovered lumi- 
nosities of artificial point-sources added to the image, 
we constructed a new "fake" response matrix. Although 
this procedure assumes that the "real" matrix is formally 
correct, the differences in the faked matrices should ac- 
curately indicate the magnitude of the statistical errors 
in that matrix. We re-fitted the composite XLF with 
appropriately added "faked" responses for each Monte- 
Carlo simulation, and assessed the l-cr scatter in the best- 
fitting parameters, indicating the level of this source of 
error. We found errors of ±0.02, 0.06x lO'^^erg s~^ and 
0.02 for Pi, Lbreak and P2, respectively, which are con- 
siderably smaller than the statistical errors reported in 
Table |. 

If the XLFs of all the galaxies in the sample are es- 
sentially the same, and it is simply the lack of counting 
statistics which have enabled us to fit single power law 
models to each dataset individually (§ ^^), we would 
expect that the best-fit composite XLF model should 
alsoorovide a good fit for each galaxy. We list in Ta- 
ble H the quality of fit for this model when applied to 
each galaxy in turn, and the best-fit models are shown 



in Figs ^-g. On a case-by-case basis, we find some sys- 
tems which marginally seem inconsistent with having a 
break at - 2 x lO^^erg s"^; for example NGC 1399. In 
contrast, the broken powerlaw model seems marginally 
preferred in some other systems, such as NGC 4552 and 
NGC 1404. However, in none of these cases is the data of 
sufficient quality to distinguish strongly between the two 
models. On average the powerlaw and broken powerlaw 
fits seem to fit the data comparably well. Nonetheless, 
if either model is 'correct', the prob(Ho) values listed for 
that model in Table || should be distributed uniformly be- 
tween and 1. We therefore tested this hypothesis with 
a K-S test. For the simple powerlaw model, the proba- 
bility that prob(Ho) are distributed uniformly is 1%. In 
contrast, for the broken powerlaw model, this probability 
is a more acceptable 65%. 

Throughout the rest of the paper, unless otherwise 
stated, we assume that the XLF of all galaxies is the same 
as our best-fitting model to the composite XLF, and ex- 
tends down to 10'^^ erg s~^. In practice, this means that 
all the 'completeness-corrected' results presented in § ^ 
refer only to LMXBs with Lx> 10^"^ erg s~^. Since the 
very low Lx XLF is dominated by only two galaxies, 
which could bias its shape if they are not representa- 
tive, we also assess how our conclusions would be af- 
fected by adopting an XLF shape which is much steeper 
in that luminosity regime (specifically an unbroken pow- 
erlaw model with (3 — 2.0). 

5. GLOBULAR CLUSTER ASSOCIATIONS 

We discuss in Paper II our data-reduction and analysis 
of archival HST WFPC2 data for 19 of the galaxies in 
our sample. In that paper, we provide detailed lists of 
all globular cluster candidates we detected, and relevant 
photometry and various derived quantities, such as GC 
specific luminosity. A number of the galaxies do not have 
WFPC2 data available, and so, for consistency, we do 
not consider them here. We consider only WFPC2 data 
for simplicity of our analysis, despite the availablity of 
superior ACS data for some of the galaxies. In practice, 
useful ACS data are available only for less than half of the 
sample, whereas most have WFPC2 data. Even with the 
relatively shallow WFPC2 data, we were able to detect, 
in most objects, more than 75% of the GC light, which 
is sufficient for our present purposes. We focus only on 
central-pointings with the WFPC2, although for some 
of the galaxies there are multiple HST pointings. We 
adopt this procedure since the Chandra PSF rapidly de- 
grades off-axis, making matching GCs and LMXBs more 
challenging, and possibly introducing a systematic bias 
between those objects with multiple pointings and those 
with only one. 

Since the absolute astrometry of Chandra is not ex- 
pected to be accurate at sub-arcsecond precision, in or- 
der to assess possible matches between LMXBs and GCs, 
we allowed a translation and a rotation transformation 
between the GC and LMXB source lists, which were ad- 
justed to maximize the number of matches. For self- 
consistency, we assumed that any GC candidate lying 
within 0.5" of a Chandra source centroid (comparable 
to the on-axis PSF) was associated with it. If, based on 
this criterion, more than one GC matches an LMXB, or if 
more than one LMXB matches a GC, we match the pair 
with the closest centroids, although this only affects a 




Fig. 8.— Composite XLFs of GC LMXBs (black) and field 
LMXBs (red), fitted with the best-fitting model for the entire 
LMXB population. 

small fraction of sources. Where both the optical and X- 
ray images of the galaxy had clearly-defined centroids, we 
constrained the transformation to ensure approximately 
the same matching precision between them as between 
GCs and LMXBs. This typically had little effect in galax- 
ies with a large number of LMXB-GC matches, but was 
an important constraint in other systems. In order to 
assess the number of spurious matches obtained by this 
procedure, we randomly generated fake GC source lists, 
scattering them approximately as the optical light until 
the total number of objects in the WFPC2 FOV matched 
that which was observed. Using the procedures outlined 
above, we determined the number of these sources ap- 
parently associated with an LMXB, which was typically 
less than ~1. 

5.1. Composite spectra of GC sources 

We next investigated potential spectral differences be- 
tween LMXBs hosted by red and by blue GCs. Exactly 
analogous to the analysis in § 3.3, we obtained compos- 



ite spectra for those LMXBs associated with red GCs 
and those associated with blue GCs. We defined any 
candidate GC with V-I>1.1 as red, and those GCs with 
V-I<1.1 as blue'^. The spectra are shown in Fig || (along 
with the composite spectrum of all LMXBs in the galax- 
ies). Both spectra were well-fitted by the best- fit models 
for the entire population (for the powerlaw models, we 
obtained xVdof=109/137 and 100/99 for the red and 
blue GCs, respectively. For the bremsstrahlung mod- 
els x^/dof was 111/137 and 102/99). There was no 
evidence of a systematically harder spectrum or intrin- 
sic absorptio n for the blue-GC source s, in conflict with 



the model of Maccaronc et al, ( |2004 ). Fitting simulta- 



neously the blue-GC and red-GC data with the same 
bremsstrahlung model but allowing A'h to vary freely 
we obtained A'h< 3 x 10^°cm"^ for the blue-GC sources 
and < 10^°cm~^ for the red-GC sources. The larger up- 
per limit for the blue-GC case in part reflects the poorer 
counting statistics but in both cases we can rule out in- 

^ The colour at which between red and blue GC sub-populations 
are separated varies from galaxy to galaxy, but we suspect the GC 
properties of LMXB hosts depend on the metallicity rather than 
the sub-population to which they belong, hence we use a fixed V-I 
division for all galaxies 
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LMXB-GC PROPERTIES 
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Note. — Details of the GC populations within the WFPC2 field of view. We show the 
total V-band luminosity in the field of view the number of GCs {Ngc), divided 

into red and blue GCs (except for NGC 1387, for which we do not have colour information; 
Paper II), the total GC luminosity {Lqc^), corrected for incompleteness, the GC source 
detection completeness, /f^, which is the fraction of the total GC luminosity detected. 
We also list the number of LMXBs observed in the WFPC2 field of view {Nx,obs), the 
(incompleteness-corrected) number of such sources in the WFPC2 field of view {Nx), the 
fraction of the LMXBs which are detected (/i^; this is derived from our fits, not strictly 
from taking Nx,obs/Nx, which is more affected by statistical noise), the number of LMXBs 
with corresponding GC counterparts {Nac,x) and the number of probable false GC-LMXB 
matches [Nfaise)- 
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Derived LMXB-GC properties 
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5.2. GC and non-GC XLF 

We accumulated composite XLFs of those LMXBs as- 
sociated with GC and those which are in the WFPC2 
fiel d bu t not identified with a GC, exactly as described 
in § 4.3. fn practice, we omitted any sources found in the 



very central regions of most of the galaxies (as discussed 
in Paper II), since features in the optical images, such 
as dust lanes, can give rise to false or ambiguous GC de- 
tections. In practice, this involved excluding regions of 
radius between 2" and 8". The XLFs were remarkably 
similar in shape and both closely resembled the compos- 
ite XLF for the entire source population (which was com- 
puted in a slightly different region). We therefore tried 
fitting both data -set s with the preferred broken powerlaw 
model fitted in § 4.3. Since the likelihood of an interloper 
being spuriously matched to a GC candidate is low, we 
assumed that all of the expected interloper sources only 
contributed to the non-GC XLF. The broken powerlaw 
model fitted both data-sets adequately (Table ||), indi- 
cating that the XLF of sources in GCs and those in the 
field are essentially the same. 

5.3. Spatial distribution of GCs and LMXBs 



We showed in § 3.4 that the spatial distribution of 



LMXBs closely follows that of red GCs (and also resem- 
bles the optical light outside the innermost ~ 0.5') in a 
subset of galaxies in the sample. If some of the LMXBs 
observed in the field are born in situ, there may, never- 
theless, be detectable differences between the radial dis- 
tributions of field and GC LMXBs. Ideally we would like 
to address whether the spatial distributions of LMXBs 
in GCs resembles that of the GC population, while for 
field LMXBs it more closely resembles the optical light. 
However, given the limited radial range over which we 
have measured the GC spatial distribution, and the rel- 
atively small numbers of sources (both GCs and, in par- 
ticular, LMXBs) that are detected, this is impractical. 
Therefore, we instead first tested, on a galaxy-by-galaxy 
basis and within the WFPC2 field of view, whether the 
radial distances of the GC and field LMXBs from the 
galaxy centre are consistent with being drawn from the 
same distribution. We assessed this with a Kolmogorov- 
Smirnov test and only considered galaxies with at least 3 
GC-LMXB matches. In all cases the hypothesis that the 
two samples were the same could not be rejected at the 
5% significance level or better. Next, we applied the same 
test to the entire population of LMXBs in the WFPC2 
field, and similarly found a null hypothesis probability of 
^--^64%, indicating excellent agreement between the two 
distributions. 

5.4. Probability that a GC contains an LMXB 

The probability that a GC is associated with an LMXB 
does not depend on the LMXB luminosity, as indicated 
by the similarities in the XLF of the GC and non-GC 
X-ray sources. Although in Paper II we compute the 
mass and metallicity of each GC, in order to investigate 
how the presence of an LMXB depends upon the GC 
properties, we focus only on the luminosity and colour. 
We adopt this approach for two reasons. Firstly, the 
metallicity errors are typically fairly large, and a non- 
negligible fraction of them have metallicities which peg 




Fig. 9. — Composite globular cluster luminosity functions of 
all GCs associated with an LMXB. The best-fit Gaussian model 
is shown (solid line), along with the best-fit model to the GCLF 
of the entire population, rescaled for clarity (dashed line). The 
models have been corrected for source detection incompleteness. 

at the maximum or minimum values allowed in our con- 
version. Secondly, as we show in Paper II, the colour 
and luminosity of the GC populations are independent 
of each other, making this analysis simpler. Since our 
optical data comprised, at most, 2 colours, we were not 
able to break the age-metallicity degeneracy and so we 
assumed that all colour depende ncv reflects a metallic- 
ity effect (c.f. Kundu et al. 2003 ). For consistency with 
past work, we make the conversion between colour and 
metallicity by using the relation of Kundu & Whitmore 
( |l998l ), log{ZFe,Gc) = -5.89 -I- 4.72(1/ - 7). To investi- 
gate now the probability that a GC contains an active 
LMXB, P(GC,X), depends upon the properties of the 
GC, we adopted the hypothesis: 



PiGC,X)^PGc~x 



^GC 



(2) 



where Lqc is the V-band luminosity of a GC, Zqc is 
its metallicity, Lq is the luminosity of the turnover in 
the globular cluster luminosity function (GCLF) for the 
entire GC population 7.2 x 10^ L©), and Zq is the 
metallicity corresponding to the peak of the GC colour 
distribution (~O.O8Z0; Paper II). The indices a and 7, 
and the normalizing constant, pgc-x, are to be deter- 
mined by the data. 

5.4.1. Globular Cluster luminosity functions 

Considering only those GC hosting an LMXB, we com- 
puted a globular cluster luminosity function (GCLF), to 
compare with that of the entire population (Paper II). 
We discuss in Paper II how the composite V-band GCLF 
is accumulated and how we fitted the data to account for 
source detection incompleteness. We assumed that the 
detection incompleteness of GCs hosting LMXBs in any 
galaxy depends upon luminosity in the same way as that 
of the whole GC population. For simplicity, we omit- 
ted NGC 1399 and NGC 1404 from this calculation, since 
they did not have any V-band photometry. We fitted the 
composite GCLF using a maximum-likelihood fitting al- 
gorithm and minimizing the Cash-C statistic exactly as 
described in Paper II. We found the GCLF of the X-ray 
luminous source population could be fitted by a single 




Fig. 10. — Composite colour distribution liistogram of all GCs 
associated with an LMXB. Also shown are the best- fit Gaussian 
model (solid line; red) and the best-fit Gaussian model for the 
entire GC population, suitably renormalized (dashed line; blue). 

Gaussian distribution, i.e. 

dnoc Ngc f {V- Vrf 
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where nac is the number density of GCs, as a func- 
tion of V-band luminosity (V), Nqc is the total num- 
ber of GCs and Vt is the apparent V-band peak lumi- 
nosity ("GCLF turnover"). Based on 100 Monte-Carlo 
simulations, we estimated prob(Ho)=46%. We found 
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in agreement with that of the entire pop- 
ulation (~1.5), but the absolute magnitude of the V- 
band turnover was -9.0to'2i approximately 2.0 magni- 
tudes brighter than th e global GCLF, c onsistent with 
previous observations ( Kim et al. 2006b ). Fig || shows 
the GCLF, along with the best-fitting Gaussian distribu- 
tion and, for comparison purposes, the best-fit Gaussian 
to the entire GC population, suitably rescaled. 

Since the re is no correlation between Lqc E^nd Zqc 
(Paper II; § 5.4.2), the dependence of P(GC,X) on these 
quantities can be assessed independently by considering 
separately the GCLF and the GC colour distributions. It 
follows from Eq || that the GCLF of X-ray hosting GCs 
will have the same shape as the entire population, but 
be systematically shifted by —0.%2aa'^ magnitudes. This 
implies that a — 1.01 ± 0.19. 

5.4.2. Globular Cluster colour distributions 

We constructed histograms of the numbers of GC as 
a function of V-I, each bin having a width of 0.05 mag. 
The histogram could be adequately fitted {via a Cash-C 
minimization algorithm) with a single Gaussian model, 
i.e. 
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where nac is the number density of GCs, as a func- 
tion of V-I colour, Ngc is the total number of GCs and 
{V — I)o is the mean of the colour distribution. We found 
{V-I)o = 1.08±0.02 mag and = 0.14±0.2 mag. This 
distribution is significantly redder and narrower than the 
GC population as a whole (Fig |l0|; Paper II), implying 
the ratio of the probability that a red GC hosts an LMXB 




Fig. 11. — Incompleteness, colour and lurmnosity corrected frac- 
tion of GCs which contain an LMXB (see EqQ), shown as a function 
of Sl. We also show the best- fitting model. 

to the probability that a blue GC does is ~2. We note 
that fitting these data with a simple Gaussian model is 
overly-simplistic since early-type galaxies generally ex- 
hibit bimodal GC colour distributions. However, the 
peaks of the metal-rich and metal-poor subpopulations 
are not cleanly separated in V-I space, and so the aver- 
age colour distributions of the clusters in a galaxy can be 
parameterized by a simple Gaussian, which is sufficiently 
accurate for our present purposes. 

The difference in a between the population of GCs as 
a whole and the subpopulation hosting LMXBs is easily 
understood in terms of the statistical errors on the pho- 
tometry, since luminous GCs preferentially host LMXBs. 
If we considered only those GCs with absolute V-band 
magnitude <-9.0, we found a — 0.14, in good agreement 
with the colour distributions of the LMXB hosts. The 
fainter GCs will have increasing statistical errors on their 
photometry, which will increase a. Assuming no corre- 
lation between the colour and magnitude of a GC (Pa- 
per II), it is trivial to show from Eq || and our adopted 
colour-metallicity relation that the colour distribution of 
GCs hosting LMXBs should have the same shape as that 
of the entire population, but be shifted by 10.87(t27 mag- 
nitudes, implying 7 = 0.33 ± 0.11. 

5.4.3. Probability that a GC hosts an LMXB 

Given the shapes of the GCLFs and GC colour distri- 
butions for our sample (for our purposes the latter can 
be approximated as a Gaussian with a peak shifted by 
Ac magnitudes from that of the entire sample, and with 
width (Tc), it follows from Eq that 



Ngc, 



X 



N 



false = PaC-X 



tTOT 

^GC j-X rGC 
T J i J i 
^0 



X exp(59.07Vc 



10.97Ac) 



(5) 



where Ngc,x is the number of GC-LMXB matches, 
N false the number of expected false matches, L^'^ the 
total luminosity of GCs, Lq the peak luminosity of the 
composite GCLF of the entire sample, is the fraction 
of LMXBs which are detected and f^'-^ is the fraction 
of the GC luminosity which is detected. For simplicity 
we here assumed that a = 1, with which our data are 
consistent. For NGC 1387 we had no colour information 
for the GCs, and so we adopted Ac = 0.0 and ctc — 0.14. 
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Fig. 12. — V-band specific frequency of LMXBs (Sn x) versus 
specific luminosity of GCs (Sl) for the WFPC2 FOV. The best-fit 
straight-Une relation is shown (solid line), along with the best-fit 
line if the intercept is fixed at (dashed line; red) 

Using the information in Tables ^ and ^, for any given 
value of pgc-x Eq || gives the number of expected 
GC-LMXB matches. Matching this model to the ob- 
served numbers of sources via the Cash-C statistic, and 
adopting a maximum-lik elih ood fitting procedure simi- 
lar to that outlined in § 4.2 we were able to place con- 
straints on pgc-x- We obtained a good fit to the data 
(prob(Hn)=63%) with pGC-x=0.025 ± 0.004. Integrat- 
ing Eq over the mean colour and luminosity distri- 
butions for the GC populations, this corresponds to a 
probability of 6.5% that a randomly chosen GC con- 
tains an LMXB with Lx> 10^^ erg s^^. We show in 
Fig 11 pgc-x computed directly from Eqn ^ for each 
dataset and the best-fit model value, as a function of GC 
specific luminosity, Sl- For comparison, if we replaced 
with the equivalent value appropriate for a power- 
law XLF, assuming j3 = 2.0, we also obtained an ac- 
ceptable fit (prob(Ho)=29%) with pGC-x=0.054 ± 0.009 
(implying 14% of GC contain an LMXB with Lx> 10^^ 
erg s^^). 

5.5. Are any LMXBs formed in the field? 

In order to test whether any LMXBs are formed in the 
field, we first adopted the simple hypothesis that: 

/ Lv 
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where Nx is the total number of LMXBs (corrected 
for source detection incompleteness), Lqq^ is the to- 
tal luminosity of the GC and Ly is the V-band lumi- 
nosity. We performed this comparison only within the 
WFPC2 FOV. For most of the galaxies in our sample, 
we tabulate in Paper II the specific luminosity of globu- 
lar clusters, defined as Sl—WOLgc/Lv- From Eq ||, ^ 
and V can be trivially obtained from a plot of Sl ver- 
sus the specific frequency of LMXBs, Sn,x, defined as 
Sn,x=8.55 X lO"^ Nx / {Ly / Lvq) (assu ming the absolut e 



V-band mag nitude of the Sun is 4.83: |Maraston| |1998D 
This is similar to the method adopted by [rwin| ( |2(305 j^ 
Adopting our preferred broken powerlaw XLF to com- 
pute Nx (Table ||), we show in Fig ^ Sn,x versus 
Sl- Taking account of errors on both x and y axes 



Fig. 13. — Incompleteness-corrected fraction of LMXBs asso- 
ciated with a GC, as a function of Sl- Shown are the preferred 
relations for the case where all sources form in GCs (solid line; 
red) and where LMXB formation in the field is proportional to the 
V-band luminosity (dashed line; blue). 



(1992), we fitted a straight line through the data. Ini- 
tially fixing 1^=0-0, i.e. requiring that all LMXBs form in 
GCs, we found a reasonable fit, (;^^/dof= 23.1/18) with 
fi = 1.20 ± 0.12- Allowing to fit freely, we found that 
the fit did not significantly improve (x^/dof=22-7/17) 
and we constrained fi = 1-271q'22 ^'^'i is < 1-8- 

For a typical galaxy with Sl~ 0-2, this implies that at 
most only ~8% of the LMXBs could be formed in the 
field- Since NGC 5845 is something of an outlier, being 
the smallest galaxy in our sample, we have investigated 
whether the results of our fit are biased by including this 
galaxy Excluding this galaxy, we obtained similar re- 
sults (xVdof=19.4/16, n = l-18l^;2i and u is < 3-4), 
implying it does not bias the fit. For comparison pur- 
poses, if we instead adopted a simple powerlaw XLF, 
with /? = 2.0, we obtained a relatively poor fit for the case 
ly = 0.0 (xVdof=39/18) with fi = 2.80t;^:^8- Allowing v 
to vary freely, the fit was similarly poor (x^/dof=38/17) 
and we constrained /i = 3.07lQ'5g and v < 3.6. 

5.6. GC LMXB retention fraction 

The probability that an X-ray binary born in a GC 
is retained by the GC, correcting for incompleteness, 
P(X,GC), is given by the relation: 
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We estimated both P(X,GC) and v by matching this 
model to the observed distribution of Ngc. x via a fit- 
ting procedure similar to that described in § ^.4.3| . Ini- 
tially we tested the hypothesis v = 0.0, i.e. all LMXBs 
form in GCs. We found that the data were well-fitted 
(prob(Ho)=81%), and we constrained P(X,GC)=0.40 ± 
0.06. Allowing v to vary freely also gave a good 
fit (prob(Ho)=95%) with V {X,GC)^Q and v = 



+ 1.7 
-5.4' 



in agreement with our results in § 5.5. Fig h 



with a strategy similar to that outlined in Press et al 
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shows values of P(X,GC) computed directly from the 
data (assuming all sources form in GCs), and the best- 
fit models- For comparison, if we adopted values of 
Nx and f^ appropriate for a single powerlaw XLF 
(with 13 = 2-0), we constrained P(X,GC)=0-55l55:i4 and 
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Fig. 14. — A comparison of our preferred cumulative XLF model 
(circles; black. The error-bars denote the approximate l-cr un- 
certainty region) with other XLFs from the literature. We have 
normalized the model to match the observed number of sources 
with Lx> ^'breafc in the Milky Way. We |mmpiirp it lyith th 



cumulative XLF of Galactic LMXBs from primm et al, 
(dashed line; red) 



(b002|) 

model tor the ditter- 



We have corrected this 
ence in energy-bands used between these authors and our work, 
under the assumption of a powerlaw spectrum ^iHth r — 9 n Jiij 



Gilfanov 



(2004 



also compare to the "universal" LMXB XLF from 

blue; dash-dot-dot-dot line), with the blue dotted lines dclmeatmg 
the approximate l-cr uncertainty region for the moHels We also 
show the best-fitting composite XLF for Kim & Fabbiano ( 2004 



light blue shaded region). The latter two models have been scaled 
rhitrarily for rlarity. We truncate both our model and that of 



|Kim fc Fabbiano above 2 X 10"^^ erg s~^. Model curves are only 
shown corresponding to luminosity ranges over which there were 
data to constrain the model. There is, in general, good agreement 
between the different models. 
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6. DISCUSSION 

6.1. The X-ray luminosity function 

6.1.1. The shape of the XLF 

When corrected for source detection incompleteness, 
we found that the XLF of early-type galaxies is very uni- 
form in shape from object to object. We found that a 
simple powerlaw fit to the XLF of each galaxy was for- 
mally acceptable, but the slope of the model correlated 
with the source detection incompleteness, implying that 
the slope flattens at low Lx- Fitting the composite XLF 
we found that it is best modeled as a broken powerlaw, 
with a break at (2.21^Q gg) x lO^^erg s~^, and negative 

differential logarithmic slopes 1.40^q '^g below the break 
and 2.84t°:3^ above. Still, a single powerlaw model with 
slope, P ~ 2.0 was able to fit the XLF of each individual 
galaxy formally ab out as well at thi s mod el. A similar 
result was found by Kim & Fabbiano (2004) for a smaller 
sample of galaxies. 

In Fig n3 we compare our best-fitting XLF, in cumu- 
lative form, with XLFs reported in the literature by a 
number of different authors. Considering the ranges for 
which each XLF has been found to be valid, we find 
excellent agreement between our best-fitting functional 
form and those reported in the lit erature, even t hat re - 



ported for Milky Way LMXBs by [Grimm et al.| ( |2002D 



We note that the absol ute normalizat ion of any of these 
models (except that of Grimm ct al.) is arbitrary. De- 



FlG. 15. — A comparison between the X-ray luminosity of point- 
sources determined from spectral-fitting (Lx^spoc) a^nd that deter- 
mined from fitting the XLF (Lx xlf) within the same region (see 
text). The best-fitting linear relation is shown. 



spite this good agreement it is still possible that the 
XLF we measure is biased at low luminosities, where the 
data were dominated by only two galaxies, NGC3115 
and NGC4552, which may not be completely represen- 
tative. Unfortunately, this cannot be ruled out with our 
current data, but will require deep XLFs measured for 
other early-type galaxies to assess it. Nonetheless, most 
of our results which depend on the XLF shape (for in- 
completeness correction) do not appear to change qual- 
itatively if we adopt an alternative XLF shape which is 
much steeper at low luminosities. 

There is some evidence in our results that there is real 
scatter in the shapes of the XLF. In order to investigate 
this further, we considered the yerv deep XLFs deter- 
mined in a few individual cases. Voss fc Gilfanov ( 2006 ) 
report an XLF for Gen A which shows a clear flatten- 



ing below 5 X 10 erg s~ , the low-Lx P is comparable 
to our measurements. Although it is difficult to inter- 
pret a cumulative XLF due to the strong correlations be- 
tween data-points, it is interesting to note that NGC 3379 
shows a significan tly flatter XLF at low Lx than (3 = 2 
( Kim ct al. 2006a), which would also be consistent with 
our broken powerlaw model. Possible counter-examples, 
however, include M 87, in which the XLF can be fltted a s 
a broken powerlaw but with pi ^ 2 (Jordan et al. 2004). 



Fitting a differential (rather than cumulative) XLi' com- 
piled from these authors' published source luminosities, 
we confirm this result, and find that our best fit to the 



composite XLF is marginally rejected. Kim et al. ( 2006a ) 
find an XLF for NGG 4278 which similarly appears to re- 
semble an unbroken powerlaw, with P ^ 2. These results 
provide further evidence that there is genuine scatter in 
the shapes of the XLFs between individual galaxies, in 
particular at low luminosities. On average, though, we 
have shown that it must flatten at low Lx- Glearly fur- 
ther, deep studies of the LMXB XLF in more galaxies 
are needed to investigate the origin of the scatter. 



6.1.2. The total luminosity of LMXBs 

In the present work, we chose to estimate the luminos- 
ity of the total LMXB populations by extrapolation down 
to 10'^^ erg s~^. This is, in part, justified by the excellent 
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Fig. 16. — Left: The integrated luminosity of the LMXBs in the galaxies, shown as a function of K-band luminosity, adopting our 
preferred, broken powerlaw XLF. We also show the mean straight line relation. Right: The deviation of Lx from its the mean relation 
shown in the left panel versus Sl- 

from the above procedure (Lx.spoc)- In a significant frac- 
tion of the galaxies, we found that Lx.spcc was extremely 
sensitive to the modelling or treatment of the background 
(consisten t with the assessment of the systematic error 
budget in Humphrey fc Buot^ 2006| ), and so these galax- 
ies have been omitted from the comparison. In our anal- 
ysis, the background is determined through a modelling 
procedure, which we have determined to be more reliable 
than adopting the background "templates" . We identi- 
fied objects in which the modelling is likely suspect, at 
least in the crucial high-energy range, by examining the 
total flux of one of the modelled spectral components — 
the powerlaw, which accounts for the unresolved frac- 
tion of the X-ray background. Although this flux is 
subject to the effects of cosmic variance, we should be 
able to resolve out a significant fraction of the brightest 
background sources (which contribute most to the cos- 
mic variance). Therefore, we should reasonably expect 
that the flux of this modelled component should be lower 
than the total flux of the AGN contribution to the back- 
ground found by, for example, Tozzi et al, (|200l|), even 
accounting for cosmic variance. 



agreement, at least to this luminosity, between our best- 
fitting composite XLF and that of LMXBs in the Milky 
Way (Fig . Another way to test the validity of this 
assumpt ion is t o ado pt a procedure similar to that out- 
lined in prwin ( 2005| ) . The ensemble average of LMXBs 
which are too laint to have been detected individually 
should contribute a hard component to the spectrum of 
the diffuse emission in each galaxy. Therefore we ex- 
tracted, for each galaxy, a spectrum from the region in 
which the XLF is computed (excluding regions around all 
detected point sources), and added the luminosity of the 
hard (unresolved LMXBs) component to the the com- 
bined luminosity of the detected sources. The details of 
the spectral extraction proc edure and the treatment o f 
the background are given in Humphrey & Buotc (2006). 

In practice, however, the procedure is rather complex 
for a number of reasons. First, some fraction of the point- 
sources are actually interloper AGNs. All sources with 
apparent Lx> 2 x 10^^ erg s~^ were assumed to be in- 
terlopers but, to minimize the effects of cosmic variance 
(which is dominated by the rare, bright point-sources), 
we disregarded these entirely. We estimated the expected 
apparent luminosity of inte rlopers over the flu x range of 
interest from the model of Tozzi et "aT] ( 2001 ) , and took 
into account the fraction of the expected flux which had 
been accounted for by the cosmic X-ray background com- 
ponent in our background modelling. Since the LMXB 
emission arises from a finite number of discrete sources, 
Poissonian statistics actually introduces an intrinsic un- 
certainty into the combined luminosity of these objects, 
irrespective of how accurately we can measure the to- 
tal Lx. We estimated the magnitude of this uncertainty 
via Monte-Carlo simulations, in which a fixed number of 
sources were randomly drawn from a luminosity distri- 
bution matching the measured XLF, and the total lumi- 
nosity measured. For each set of simulations, we were 
able to estimate the scatter (standard deviation) in the 
total luminosity as a function of its expected value. For 
each galaxy, the scatter was added in quadrature with 
the measurement errors. 

Fig |l5| shows a comparison of the luminosity deter- 
mined irom XLF fitting (Lx.xlf) and that determined 



In half of the galaxies, 
we actually found that its flux was considerably higher 
than expected from this argument, typically indicating 
systems in which our models do not completely capture 
the shape of the high-energy background. In such cir- 
cumstances, Lx,spcc is almost certainly in error (even 
though the luminosity of the hot gas emission, which has 
a much softer spectrum, is fairly well-determined), and 
so is omitted here. 

We fitted a model of the form Lx.xlf — CLx,spcc to 
these data. Since the XLF certainly extends to fainter 
LMXBs than lO^^erg s~^, we do not expect C to be ex- 
actly 1; its precise value depends on both the shape of 
the XLF and the actual low-Lx cut for LMXBs. Extrap- 
olating our preferred XLF, for example, we expect ^ to 
be in the range 0.9-1.0. In contrast, for pure powerlaw 
XLF with P = 2.0, C may be as low as 0.70 if the XLF is 
unbroken down to 10'^^ erg s~^. Computing Lx,xlf for 
each galaxy with our best-fitting composite XLF, we ob- 
tained a good fit with this model (x^/dof=l 1.1/11), with 
C = 0.87 ± 0.12, which is consistent with expectations. 



Adopting Lx.xLF from the powerlaw fits to the XLF, we 
found C — 1.06 ± 0.16, which impHes that the XLF must 
be truncated below a Hmit which is >5 x 10'^^ erg s~^. 

Previous authors have reported a strong correlation 
between the total luminosity of the X-ray sources and the 



K-ba .nd lumino sitv of the galaxy (Gilfanov 2004; Kim & 
Fabtiano 2004). Using our larger sample or galaxies, we 
find a similar correlation (Fig |l6| ), corresponding to 

LLMXB/Lj, ^ 0.11 X lO^Oergs-iL-^ (8) 

with a scatter of 0.06 x lO'^'^ergs^^Lj^^, assuming our 
preferred XLF shape for each galaxy. For comparison, if 
we assumed that the XLF is a powerlaw with (3 = 2.0, 
the mean and scatter (omitting IC4296, for which the 
error-bars are large) became 0.17 x 10'^*'ergs~^Lj^Q and 



0.09 X 10 ergs L^^, respe c tively , in excellent agree- 
ment with Kim fc Fabbiano (2004). In both cases the 
statistical errors are considerably smaller than the in- 
trinsic scatter. 

At face value this correlation indicates that the inte- 
grated LMXB populatio n is a good ind icator of the stel- 
lar mass i n the galaxy ( Gilfanov 2004). However, as we 
show in § |3.5| , the majority, if not all, of the LMXBs 
must be formed in GCs. The relation between the num- 
bers of LMXBs and the combined luminosity of the GCs 
does not show any scatter larger than the (albeit large) 
statistical errors. There is, however, a correlation be- 
tween luminosity and the numbers of GCs in a galaxy 
suggesting that the correlation found in Fig |l^ may arise 
through this third parameter. We show in the right panel 
of Fig 16 the deviation of Lx from the mean relation 
with respect to Lk as a function of Sl. This shows a 
clear correlation, which would support this interpreta- 
tion, although this figure should be treated with caution 
since Sl and Lx are not computed from self-consistent 
apertures. 

6.1.3. The high-Ijyi slope 

Given the quality of the data, it was not possible to in- 
vestigate how Lbreak Or f3i vary between galaxies. How- 
ever, we were able to investigate whether P2, the high-Lx 
slope, does so, possibly providing further clues as to the 
origin of the XLF shape. Where it could be constrained, 
we allowed (32 to fit freely and the results are shown in 
Table We searched for correlations between it and var- 
ious interesting parameters such as age, metallicity and 
[a/Fe] for the stellar population, Sl and the peak colour 
of the GC distribution. The ages and met allicitics listed 
in Tabl e [l| w ere c omputed as desc ribed in Humphrey & 
Buote (2006) and Humphrey et al. ( 2006| ). For those sys- 
tems for which the age had been fixed to 12 Gyr in our 
previous analysis we relaxed this constraint when esti- 
mating the age. 

If the high-luminosity sources are domi nated bv s hort- 
lived objects forming in the field (e.g. Wu 2001), we 
would expec t to see a correlation betw een /32 and the 
stellar age. Ivanova & Kalogera (2006) constructed a 
model for the high-Lx XLF slope, which they related to 
the transient duty cycle. In this picture, if mass-transfer 
is dominated by red giant or white dwarf donors, P2 
should also depend on age. However, we did not find 
any evidence of a correlation (prob(no correlation)=70% 
for Kendal's r-test, and we obtained similar results with 
Spearman's rank test) between P2 and age (Fig |l7|). 



In contrast we did find a marginally significant corre- 
lation between [a/Fe] and P2, as shown in Fig |^ for 
which prob(no correlation)=3-8% for the different tests. 
To investigate whether such a correlation is robust in 
the presence of the large statistical errors, we investi- 
gated the effect of randomly scattering the data-points 
within their error-bars. We found that 56% of the re- 
sulting simulations gave a significant correlation (in com- 
parison to the 5% expected if the correlation is not ro- 
bust). Given the number of different trials we performed, 
the correlation between [a/Fe] and P2 should be consid- 
ered tentative at best, and needs to be confirmed with 
a larger data-set. Nonetheless, taking it at face value, 
it is intriguing. One possible cause of such a correlation 
is a systematically- varying IMF. A more top-heavy IMF 
would increase [a/Fe] and also affect the compact object 
demographics. If there are more massive objects, it is 
unsurprising if more sources are observed which are able 
to sustain high Lx. Alternatively, enhanced [a/Fe] in the 
stars may significantly affect how mass transfer operates, 
increasing the population of X-ray bright LMXBs. How- 
ever, neither the peak GC colour nor overall metallicity 
appear to correlate with P2, and it is not clear why [a/Fe] 
should have a stronger impact on the accretion process 
than total metallicity. Clearly both of these scenarios re- 
quire that either the high-Lx sources are systematically 
biased to form in the field (since most LMXBs do not), 
or the variations in stellar [a/Fe] (and/ or the IMF) must 
also be reflected in the GC populations. 

6.1.4. The meaning of the break 

Early results suggesting the presence of an XLF break 
at ~ 2 X 10^^ erg s~^ have been called into question 
due to the failure of the authors to correct properly fo r 
source detection incompleteness (Kim & Fabbiano 2004). 
Although we now confirm the presence ol such a break 
even when the data have been corrected for incomplete- 
ness, a break was not required statistically in any indi- 
vidual galaxy. The significance of the feature in some 
early analyses may have been exaggerated, therefore, by 
incompleteness effects. 

It has been suggested ( Sarazin et alj |200l| ) that the 
break represents a division between neutron-star and 
black-hole binary systems, with objects above this limit 
being entirely black hole binaries, and those below the 
limit a mixture of black hole binaries and neutron-stars. 
It is difficult to reconcile this hypothesis, however, with 
the shape of the XLF. Most systems are not observed 
close to their Eddington limit; for example all neutron- 
star LMXBs do not have Lx-- 2 x lO^^erg s~^ There- 
fore, the black-hole binary XLF would be expected to 
extend below this limit and, crucially, would not be ex- 
pected to exhibit any feature at this characteristic lu- 
minosity. An extrapolation of the high-Lx slope (com- 
prising only the putative black hole population) dramat- 
ically over-estimates the numbers of sources at low Lx, 
since the XLF flattens below t he break, effec t ively rul- 
ing out this picture. Although Barazin et al. (2001) ac- 
tually recognized this problem, they argued there must 
be different source populations either side of the break. 
However, this would still, in fact, require the black hole 
binary XLF to "care" about the Ed dington limit of a 
neutro n-star, which is not expected. Bildsten fc Dcloye 



( 2004 ) explained both the break and low-Lx slope of the 
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Fig. 17. — The slope of the XLF at Lx> 2 x 10^® erg s"-"^ versus age (left panel) and [a]/Fe (right panel) of the stellar populations. 



XLF found by |Kim fc Fabbiano| (|200J) in terms of ultra- 
compact X-ray binaries. However, the break luminos- 
ity we detect is significantly lower and, more problem- 
atically, the low-Lx slope predicted {f3 ~ 1.8) is much 
steeper than we detect. 

It is worth noting that a further problem with inter- 
preting the break as the Eddington limit is, of course, 
that the point source Lx we measure in the Chandra band 
is certainly an underestimate of the bolometric luminos- 
ity. For example, for the representative Galactic neutron- 
star X-ray binaries X 1916-053 and X 1624-490, an ex- 
trapolation of the best- fitting models for broad-band 



Bep poSAX observa tions ( Church et al. 1998; Balucihska- 



Chui |ch et al.||2000 ) imphes the 0.3-7.0 keV flux under 
estimates that in the 0.1-200 keV band by factors ~2.8 
and 1.5, respectively, and obviously the bolometric fiux 
is underestimated by even larger factors. 

Another possible explanation for the break is that it 
arises from an evolving LMXBs popula tion born in a 



utative recent burst of star formation (Jeltema et al 



E003; Wu 2001). The absence of any obvious correlation 
between the high-Lx shape and the age of the stars seems 
to rule out this model. Furthermore, given the short 
lifetime of an X-ray binary (< lO^yr) this model seems 
unli kelv to be able to sustain a break at 2 x lO'^^erg s~ ^ 



Way have Lx~l-3xl038erg s'^ ( |Grimm et "aL| |2002| ), 
measured in the Chandra band, in good agreement with 
Lbreak in the early-type galaxy XLF. If, as in the Milky 
Way, sources brighter than this limit are predominantly 
lower-Lx sources exhibiting brief periods of outburst or 
flaring, this would naturally cause a break. This can be 
explicitly tested by looking for characteristic patterns of 
(spectral) variability in the brightest LMXBs. In the cur- 
rent work there were, unfortunately, insufficient photons 
to allow us to do this. One possibl e coun t er-arg ument to 
this picture is the observation by Irwin (2006) that no 
transient behaviour was seen in the brightest sources in 
NGC 1399 and M 87, over a timescale of ~a few years, 
effectively placing a limit on transient behaviour in these 
objects of ~50 yr. Although these objects are obviously 
not short-term transients, in no way does this work sug- 
gest they are not transient over periods of a few decades 
or longer, and so our interpretation still holds. In any 
case, such long-term transients appear very rare (as seen 
in the Milky Way). 

As a corollary, it is interesting to note that source vari- 
ability may have other impacts on the shape of the XLF. 
The characteristic patterns of short term variability (e.g. 
White et al. 1995 ) are superimposed on long term sec- 



in a typical galaxy with a ^12 Gyr stellar population. lowski et al, 2002 ) 



ular evolution in the mass-transfer rate (e.g. Podsiad 
so that the locus of an individual 



In any case, we did not find any convincing evidence 
of a break which systematically evolves from galaxy-to- 
galaxy and most, if not all, LMXBs are formed in GCs. 

The well-studied population of Milky Way LMXBs 
provides some intriguing clues as to what may be the 
origin of the break. Galactic LMXBs are known to 
be h ighly variable (e.g. |Liu et a~ |200l| ; Grimm et al. 
002). In particular, all confirmed LMXB black hole 



ina ries in the Milky Way are transient sources ( Mc 



m 

Clin qock fc Remillarq |2005| ) , and there is evidence that 
a significant portion of LM XBs in early-type galaxies 
may also be transient (e.g. Kraft et al. 2001). Differ- 



source in the XLF should evolve with time. It is pos- 
sible, therefore, that the time-averaged lightcurve may 
be the dominant factor in determining the shape of the 
XLF, which might explain the comparative uniformity 
in its shape from galaxy to galaxy. Obviously compact 
object demographics must play a role since black-hole bi- 
nary systems exhibit, on averag e, slightlv different prop - 
erties to neutron-star systems ( Fanaka &: Lewin 1995| ). 



ent classes of LMXB, which exhibit strikingly different 
patt erns of variability occupy distin ct luminosity bands 
(e.g. Hasinger & van der Klit 1989), so that the XLF is 
not expected to be "scale-free . It is striking that the 
so-called "Z-track" sources, which comprise the majority 
of the most luminous persistent LMXBs in the Milky 



Therefore differences in the IMF might still produce the 
possible correlation we see between [a/Fe] and ^2- 

6.2. The GC-LMXB connection 

We find a strong correlation between the number of 
GCs in the galaxy and the number of LMXBs, placing 
tight constraints on the fraction of sources formed in the 
field. We find that the number of LMXBs formed in 
the field is < 1.8 x 10-^°Lv/Lvq, implying < 10% of 
the LMXBs in a galaxy with 0.2. In fact, the data 



are consistent with the hypothesis that all the LMXBs 
are formed within GCs. IdeaUy we would have per- 
formed this comparison using the K-band, rather than 
the V-band hght, since it better traces the mass of the 
underlying stellar population. However, since our best- 
fitting relation passes through the origin, this indicates 
the number of LMXBs is consistent with being exactly 
proportional to the combined V-band luminosity of the 
GCs. Since the GCs are expected to be uniformly old 
stellar systems (e.g.. Paper II), the V-band mass-to- light 
ratio should not vary considerably between the galaxies. 
Scaling both the number of LMXBs and the luminosity 
of GCs for each galaxy by the same arbitrary factor (e.g. 
the K-band luminosity of the galaxy; hence computing 
Sn,x and Sl using the K-band, rather than V-band lumi- 
nosity) should not affect the quality of a "yc>cx" fit to the 
data. Although the exact numerical relations we derive 
depend on the shape of the XLF we adopt for incomplete- 
ness correction, qualitatively our conclusions are found 
to be insensitive to the exact XLF shape. 

Our c o nclus ions are significantly different from those 
of Irwin (|2005|) , who, using a smaller sample of galaxies, 
argued that the total luminosity of the LMXBs divided 
by the optical luminosity of the galaxy is proportional 
to Sn plus a constant offset. We suspect the reason for 
this discrepancy arises from the way in which t he lu- 
minos ity of undetected LMXBs was estimated. [Irwin 



( p005| ) 



computed the luminosity of unresolved LMXBs 
from the high-energy (2.0-6.0 keV) portion of the diffuse 
emission spectrum in each galaxy (excepting NGC 1399). 
This implicitly assumes that there is no contribution 
from the hot ISM at these energies, which for a num- 
ber of the galaxies in the sarnple ruay not be strictly 
true (e.g. Humphrey fc Buotc 2006 ; Humphrey et al 



^0041) 



More problematically, the determination of the 
unresolved source flux from the spectrum is very sensi- 
tive to the tr eatment of the background (§ 6.1.2 ). Al- 
though Irwin took a local background from the outer 
parts of the S3 chip, this fundamentally requires that 
there is no extended diffuse emission in any of the galax- 
ies, which is not true for all of the galaxies in the sample 



(e.g. Humphrey et al. 2004), and requires careful treat- 



ment of vignetting effects. In fact, he found that, when 
attempting spectral fitting, the models systematically 
underestimated the high-energy spectrum, in contrast 
t o our results when carefuU y modeling the background 
(Humphrey & Buote 2006), suggesting that the back- 
ground may not have been completely accounted for. 

We found that '--^l currently active (Lx> 10'^''erg s~^) 
LMXB was formed per IO^Lq V-band integrated lu- 
minosity of GCs, of which -40-50% of the LMXBs 
are retained within the GC. Approximately 7% of 
the GCs within the WFPC2 field contain an active 
LMXB, in rough agreement with the fr action seen 
in the Milky Way and, possibly, M31 ( pi Stefano 
ct a[] |2002| ). If the hfetimc of an LMXB is tl 
and the fraction of the time it is active with Lx> 
lO^^erg s~^ is Fd, our results imply an LMXB forma- 
tion rate in GCs of 



lO-i^L-J^yr 

-1 

d 



1.5Gyr-i(TL/10V)"^F 



HrL/10V)-^F, 



per GC, averaging over 



or 

the mean GCLF. 

Another piece of evidence which supports the hypoth- 
esis that all LMXBs form in GCs is the good agreement 



between the spatial distribution of the LMXBs and the 
GCs, when one considers an appropriately-weighted com- 
bination of the red and blue GCs. We find no evidence of 
a difference in the spatial distributions of LMXBs found 
in GCs and the field implying that field LMXBs are most 
likely ejected only with modest kick yelocities fr or n GCs 



With a smaller sample of galaxies Kim et al. ( ^006b ) 



found similarities in the distributions of the LMXBs in 
the field and in GCs, consistent with our result, but found 
that the LMXB distribution was far steeper than the 
GCs. These authors' optical data covered a larger radial 
range than ours, and so it is possible the agreement we 
find with our data is an artifact of our smaller field of 
view. However, we note that these authors did not cor- 
rect for the spatial variation in source detection incom- 
pleteness, which affects radial profiles of both the LMXBs 
and the GCs, and may have played a r ole in this discrep- 
ancy. In contrast Kundu et al. (2007) report, based on 
a sample of 5 galaxies, that the LMXB-hosting GCs are 
distributed similarly to the GCs, whereas field LMXBs 
have a more centrally-peaked distribution. To some ex- 
tent, we suspect that small number statistics also plays 
a role in the discrepancies between thes e auth ors' work 
and ours. The sample of Kundu et al. (2007), for ex- 
ample, contains the SO galaxy NGC 3115, which had the 
most centrally-peaked LMXB distribution in our sample. 
Furthermore, these authors compare the distribution of 
LMXBs and that of all the GCs, which is more extended 
than an appropriately-weighted combination of red and 
blue GCs (see §U). 

Having corrected for source detection incompleteness 
in both the LMXB and GC populations, we found lit- 
tle evidence that the fraction of LMXBs coincident with 
GCs vary from g alaxy to gal axy. This is in sharp contrast 
to past studies ( [luctt 2005), which argued for a correla- 
tion between this fraction and the specific frequency of 
globular clusters, which implies a significant fraction of 
the LMXBs are formed in the field. This study, how- 
ever, drew upon literature results which did not always 
treat the data consistently. In particular, literature val- 
ues of Sn were used, which were typically not computed 
in the same aperture as was used for LMXB-GC match- 
ing. Since the distribution of GCs as a whole is fiatter 
than the optical light (or the distribution of LMXBs, 
which follow the more centrally-peaked red GC distribu- 
tion), it is crucial that Sn is computed within this aper- 
ture. Adopti ng loca lly-computed specific frequencies the 
cor relation of J juett is seriously degraded, as pointed out 
by Kim et aLf ( 2006b| ). Furthermore, when one corrects 
for incompleteness, there is little evidence of any statis- 
tically significant variation from galaxy-to-galaxy. 

We found that the probability a given GC contains an 
LMXB is proportional to its luminosity to the power of 
1.01±0.19 and its metallicity to the power of 0.33±0.11. 
These values are v ery s imilar to recent results obtained 
by ISivakoff ct al| ( |2007|), [Jordan ct al.| (|00|, for M 87) 



and the sample of 3mits ct al. "( 2006^ . In order to in- 



terpret these results in terms of physical models, it is 
necessary to recast Eqn ^ in terms of Mqc, the mass of 
the GC. For the range of stellar metallicities implied by 
the colours of LMXB-hosting GCs {[Fe/H] ~-l.6-0.25), 
the stellar M/Ly ratio actually depends on metallicity. 
To illustrate this point, we adopted M/Ly ratios for an 
assumed 10 Gyr population with a Kroupa IMF from 



MarastonI ^99§)* and found M/Ly ^ 3.5Z'pl^, which 
is accurate to better than ^10% for the range of inter- 
est. Inserting this into Eqn 0, we obtained P{GC, X) ~ 
0.023(Mgc/Mo)1-"i±o-19(Zgc/^o)°-"*°-'', where Mq = 
1.6 X lO'^M^T,. 

Maccarone et "aT" (2004) proposed an irradiation- 



LMXBs in red and blue GCs have essentially the 
same spe ctrum, ruling out the ir radiation induced 



model of Maccarone et al. (2004) 



induced wind model to account for an excess of red GCs 
harbouring LMXBs, which predicts an approximately 
powerlaw dependence on metallicity, with an exponent 
^0.39. This is stronger than the dependence we obtained 
(differing at ~ 2a). Furthermore, we note that the pre- 
dicted enhanced absorption in the blue-GC LMXBs than 
red-GC LMXBs i s not observed in our dat a. It is interest- 
ing to note that Maccarone et al. ( 2004| ) also estimated 
that the effects of metallicity on stellar radius should 
lead to a dependence P{GC,X) oc Zq^, which is in 
good agreement with our estimate above. We stress that 
our value of the metallicity exponent is derived assum- 
ing that all of the colour variation between GCs arises 
solely from metallicity effects, ignoring the effects of age, 
which would affect both the M/L ratio and the conver- 
sion from GC colour to metallicity. Clearly further multi- 
colour or deep spectroscopic observations to break the 
age-metallicity degeneracy are needed to investigate this 
further. 

7. SUMMARY AND CONCLUSIONS 

Using a sample of 24 galaxies observed with Chandra 
we found, in summary: 

1. Correcting for source detection incompleteness, the 
point-source XLF of individual galaxies is consis- 
tent with a single powerlaw, with f3 ~2.0. We find 
that P correlates with incompleteness, indicating 
that the XLF steepens at low Lx. 

2. The composite XLF of all the galaxies is best- 
fitted by a broken powerlaw with a break at 
2.2lto fe X lO^^erg s"! and slopes 1.401^ !^] and 



-0.13 



2.84^Q 3Q above and below the break, respectively. 

3. The shape of the XLF is inconsistent with the break 
representing a strict division between black hole 
and neutron-star systems. However, in common 
with the Milky Way, it may represent the lumi- 
nosity of the brightest persistent LMXB, so that 
sources with higher Lx are exhibiting flaring or in 
transient outburst, naturally producing a break. 

4. The slope of the XLF at high luminosities does not 
correlate with the stellar population age, but shows 
evidence of a weak correlation with [a/Fe], suggest- 
ing a possible effect of the IMF on the XLF. 

5. We found that the combined LMXB luminos- 
ity is approximately proportional to Lk of the 
galaxy. However, there is significant intrin- 
sic scatter in the relation, such that L]^^^ /Lk 
has a mean and standard deviation of 0.11 and 
0.06 X 103"ergs-iLKQ. 

6. We find no difference in the XLF or composite spec- 
tra between LMXBs in GCs and those in the field. 



7. Correcting for incompleteness, GC colour and 
GCLF effects, we find that the probability a GC 
hosts an LMXB with Lx> 10^^ erg s^^ is 6.5%, 
and the probability an LMXB is in a GC is ~40%. 
These do not vary from galaxy to galaxy. 

8. The probability that a GC hosts an LMXB depends 
on luminosity and metallicity to the powers a — 
1.01±0.19and7 = 0.33 ±0.11, respectively. When 
the metallicity dependence of the stellar M/L ratio 
is taken into account, the tendency for red GCs 
preferentially to host LMXBs may be consistent 
solely with metallicity infiuences on stellar radii. 

9. The specific frequency of LMXBs is proportional 
to the Sl without a significant constant off- 
set, implying that all LMXBs form in GCs. 
Our results imply an LMXB formation rate of 
- 1.5Gyr-i(TL/10Vr)"^F-i per GC (where tl is 
the LMXB lifetime and the fraction of its life 
spent with Lx> 10'^^ erg s~^), and put an up- 
per limit of - 1.8 X IQ-'^^Lv/Lvq LMXB with 
Lx> 10"^^ erg s^-^ formed in the field. 



10. This is further supported by the fact that the 
LMXB distribution closely follows the distribution 
of red GCs, when incompleteness is taken into ac- 
count. Furthermore, the spatial distributions of 
GC and field LMXBs are essentially the same. 

11. With the present data we cannot completely rule 
out that the two galaxies which dominate the XLF 
at the lowest Lx are not representative, which 
could distort the XLF shape in that region. Al- 
though this would have an impact on our numerical 
results, qualitatively our conclusions appear insen- 
sitive to the low-Lx slope. 
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* Using the updated model-grids made available by the 
author at http://www-astro.physics.ox.ac.uk/~maraston/ Clau- 



dia's_Stellar _Population_Models.html 



APPENDIX 
SOURCE LISTS 

We include in Table || the point-source lists of all the galaxies in the sample. We include only those sources within 
D25 of each galaxy, which should mitigate against contamination by background AGN. For each source, we provide 
the source name, which incorporates the coordinate information , th e total number of counts detected, the 0.3-7.0 keV 
X-ray luminosity (Lx), the two hardness ratios HRl and HR2 (§ 3.1) and the projected distance from the galaxy centre 
(Afl). For those LMXBs within the appropriate WFPC2 field, we list under "GC" the GC catalogue number given 
in Paper II for any which match GCs. If the source has no match, we report "None" , and if the source is not in the 
WFPC2 field, it is marked by an ellipsis. 
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